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1  Introduction 


For  many  years,  conventional  radars  transmitted,  received,  and  processed  the  same 
waveform  on  every  pulse  or  burst  within  a  coherent  processing  interval  (CPI),  independently 
of  the  environment.  In  the  70s,  adaptive  processing  began  to  be  developed  and  radars  began  to 
be  more  flexible  on  receive.  For  the  first  time,  the  processing  of  received  signals  changed 
depending  on  the  environment  (noise,  clutter,  interferences). 

Now,  modern  radar  systems  have  considerable  flexibility  in  their  modes  of  operation,  both  on 
receive  and  transmit.  In  particular,  it  is  possible  to  modify  the  waveform  on  a  pulse  to  pulse 
basis,  and  electronically  steered  phased  arrays  can  quickly  point  the  radar  beam  in  any 
feasible  direction.  Such  flexibility  calls  for  new  methods  of  designing  and  scheduling  the 
waveforms  to  optimize  the  radar  performance.  Then,  an  agile  and  diverse  waveform  radar 
system  should  be  able  to  change  on  the  fly  the  transmitted  waveform  based  on  the  information 
estimated  or  a  priori  known  on  the  environment,  on  the  targets  and  the  jammers. 

Moreover,  in  a  radar  network  each  sensor  should  also  be  able  to  operate  and  perform  its  task 
without  negatively  interfering  with  the  other  sensors  and,  possibly,  to  improve  the 
performance  of  the  whole  network.  Then,  the  waveforms  used  by  the  radars  in  a  complex 
network  should  be  designed  and  changed  on  the  fly,  based  on  the  clutter,  target  and 
interference  echoes;  they  should  guarantee  good  target  detection  and  parameter  estimation  in 
different  scenarios  and  should  allow  an  optimal  access  to  the  same  transmit  channel. 

In  the  first  report  of  the  research  activity  [TR1]  we  introduced  the  Ambiguity  Function,  an 
analytical  tool  for  waveform  design  and  analysis  that  is  useful  for  examining  resolution, 
sidelobe  behavior,  and  ambiguity  in  range  and  Doppler  of  a  given  signal  waveform.  As  well 
known,  to  improve  resolution  without  worsening  detection  performance,  radars  should  use 
pulse  compressed  waveforms.  These  waveforms  are  obtained  by  adding  frequency  or  phase 
modulation  to  a  rectangular  pulse.  There  are  a  vast  number  of  pulse  compressed  waveforms  in 
the  literature;  the  most  commonly  used  from  the  radar  community  are  the  Frequency  Coded 
Waveforms  (FCW),  that  has  been  analyzed  in  detail  in  [TR1].  These  kind  of  waveform  are 
often  used  in  high  resolution  radar  systems  because  they  guarantee  good  target  detection  and 
estimation.  FCW  signals  are  characterized  by  an  auto-ambiguity  function  that  exhibits  a 
narrow  thumb  tack  shape  with  low  sidelobes.  In  contrast,  in  applications  like  multi-access 
communications,  attention  is  paid  in  designing  a  sequence  of  frequency  coded  waveforms 
with  small  cross-correlation  functions.  In  multi-user  radar  system  scenarios  both  objectives 
are  desirable.  Unfortunately,  there  is  a  tradeoff  between  these  objectives.  Frequency  hop  pulse 
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trains  based  on  Costas  codes,  for  instance,  are  known  to  have  almost  ideal  auto-ambiguity 
function  but  not  very  good  cross-ambiguity  properties.  On  the  contrary,  frequency  coded 
signals  based  on  linear  congruences  have  ideal  cross-ambiguity  properties  but  unattractive 
auto-ambiguity  properties.  Some  techniques  to  design  multiple  access  frequency  hop  codes 
with  good  auto  and  cross-ambiguity  functions  have  been  described  in  [TR1];  in  our  work  we 
focused  on  the  frequency  hop  patterns  constructed  upon  an  extension  of  the  theory  of 
quadratic  congruences  and  we  analyzed  also  a  scenario  composed  of  two  radars  transmitting 
in  the  same  band,  which  can  illuminate  the  same  area  looking  for  same  target.  In  particular, 
these  systems  can  use  either  the  same  or  different  transmitted  FCW.  In  [TR1]  we  investigated 
the  impact  of  the  presence  of  the  transmitted  signal  of  the  second  radar  (the  interfering  radar) 
on  the  first  one  (the  reference  radar)  in  the  estimation  of  the  target  direction  of  arrival  (DOA). 
The  considered  target  DOA  estimators  were  the  Monopulse,  the  Pseudo  Monopulse  (PM),  and 
the  Maximum  Likelihood  (ML)  techniques.  The  performances  of  these  estimators  have  been 
numerically  evaluated  in  terms  of  root  mean  square  error  (RMSE)  and  probability  of  loosing 
the  target  in  different  operational  conditions.  Clearly,  the  performances  of  the  reference  radar 
are  strongly  dependent  on  the  codes  used  by  the  two  radars;  therefore,  using  proper  coding  the 
influence  of  the  interfering  radar  can  be  significantly  reduced.  This  was  the  main  topic  of 
[TR2]  where  we  presented  two  modified  version  of  the  PM  and  ML  estimators,  the  Modified 
Pseudo-Monopulse  (MPM),  the  Modified  Maximum  Likelihood  (MML),  the  Sub- sequent 
Pseudo-Monopulse  (SPM)  and  the  Sub-sequent  Maximum  Likelihood  (SML).  In  this  report 
we  describe  briefly  only  the  MPM  technique  for  the  lack  of  space,  for  more  details  we  refer 
the  reader  to  [TR2].  In  [TR1]  and  [TR2]  we  analyzed  a  netted  system  composed  by 
monostatic  radar  system.  In  the  most  general  case,  in  a  multistatic  radar  system  the  receivers 
can  also  exploit  the  signal  emitted  by  non  colocated  transmitters.  To  this  end,  it  is  very 
important  to  define  and  analyze  bistatic  radar  system.  This  was  the  main  topic  of  [TR3]  and 
[TR4].  Bistatic  radar  may  be  defined  as  a  radar  in  which  the  transmitter  and  receiver  are  at 
separate  locations.  In  monostatic  configuration,  estimation  of  the  time  delay  and  Doppler  shift 
directly  provides  information  on  range  and  velocity  of  the  target.  This  is  possible  also  in  the 
case  of  bistatic  radar  configuration,  even  if  the  relation  between  measured  or  estimated  time 
delay  and  Doppler  frequency,  and  target  distance  and  velocity  is  not  linear.  To  measure  the 
possible  global  resolution  and  large  error  properties  of  the  target  parameters  estimates,  the 
ambiguity  function  (AF)  is  often  used  both  in  mono  and  in  multistatic  scenarios. 
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We  also  exploit  the  relation  between  the  AF  and  the  Cramer- Rao  lower  bound  (CRLB)  to 
calculate  the  bistatic  CRLBs  of  target  range  and  velocity.  The  bistatic  CRLBs  provide  a  local 
measure  of  the  estimation  accuracy  of  these  parameters.  Moreover,  we  compare  monostatic 
and  bistatic  CRLBs  as  a  function  of  the  number  of  integrated  pulses,  target  direction  of  arrival 
(DOA),  and  bistatic  baseline  length.  The  information  gained  through  the  calculation  of  the 
bistatic  CRLBs  can  be  used  in  a  multistatic  radar  system  for  the  choice  of  the  optimum 
transmit-receive  pair.  Multistatic  radars  utilize  multiple  transmitters  and  receivers.  Such 
systems  differ  from  typical  modem  active  radars  since  they  consist  of  several  different 
monostatic  and  bistatic  channels  of  observation.  Due  to  this  spatial  diversity,  these  systems 
present  challenges  in  managing  their  operation  as  well  as  in  usefully  combining  the  data  from 
multiple  sources  of  information  on  a  particular  area  of  surveillance.  The  information  gain, 
obtained  through  this  spatial  diversity,  combined  with  some  level  of  data  fusion,  can  give  rise 
to  a  number  of  advantages  over  both  the  individual  monostatic  and  bistatic  cases  for  typical 
radar  functions,  such  as  detection,  parameter  estimation,  tracking  and  identification.  As 
showed  in  [TR3],  the  performance  of  each  channel  of  the  multistatic  system  heavily  depends 
on  the  transmitted  waveform  and  on  the  geometry  of  the  scenario,  that  is,  the  position  of 
receivers  and  transmitters  with  respect  to  the  position  of  the  target.  In  particular,  both 
geometry  factors  and  transmitted  waveforms  play  an  important  role  in  the  shape  of  the 
Ambiguity  Function  and  hence  in  the  value  of  the  CRLBs.  In  [TR4],  after  we  calculated  the 
bistatic  CRLBs  of  target  range  and  velocity  of  each  TX-RX  pair  as  a  function  of  the  target 
kinematic  parameters  and  to  provide  a  local  measure  of  the  estimation  accuracy  of  these 
parameters.  In  [TR4],  we  exploited  the  results  obtained  from  this  analysis  to  compute  the 
rules  for  selecting  the  best  weighting  coefficients  for  fusing  the  signals  from  multiple 
receivers  in  order  to  improve  the  detection  performance  and  the  estimation  accuracy  of  the 
kinematic  parameters  of  the  target.  We  also  introduced  an  optimization  methodology  for 
selecting  only  some  channels  for  the  network,  independent  of  the  adopted  fusion  rule. 

The  results  derived  in  [TR3]  and  [TR4]  can  also  be  used  to  analyzed  the  performance  of  a 
Passive  Coherent  Location  (PCL)  system,  that  is,  bistatic  radars  that  exploit  the  signal  emitted 
by  illuminators  of  opportunity,  such  as  broadcast  or  communications  signals.  PCL  systems 
have  some  significant  attractions,  in  addition  to  those  common  to  all  bistatic  radars:  the  PCL 
receivers  do  not  need  any  transmitter  hardware  of  their  own,  are  completely  passive,  and 
hence  undetectable.  In  [TR5]  we  calculated  the  monostatic  and  bistatic  Ambiguity  Function 
(AF)  of  a  Quadrature  Phase  Shift  Keying  (QPSK)  signal  where  the  pulses  are  shaped  with  a 
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Root  Raised  Cosine  (RRC)  filter.  The  monostatic  and  bistatic  modified  Cramer-Rao  lower 
bounds  (MCRLBs)  for  the  estimation  of  target  range  and  velocity  are  also  derived  and 
analyzed.  The  QPSK  modulation  is  used  in  the  downlink  of  a  Universal  Mobile 
Telecommunications  System  (UMTS)  base  station,  hence  the  results  of  our  analysis  provide  a 
useful  tool  to  asses  the  performance  of  a  Passive  Coherent  Location  (PCL)  system  where  the 
non  co-operative  transmitter  of  opportunity  is  a  UMTS-base  station.  The  actual  growing 
coverage  of  UMTS  signals  on  the  international  territory  makes  multistatic  radar  configuration 
feasible,  therefore  these  results  can  also  be  exploited  for  the  dynamical  selection  of  the 
transmitter  in  a  multistatic  radar  system  where  multiple  UMTS  base  stations  are  used. 

All  the  obtained  results  are  briefly  summarized  in  this  report,  for  more  details  we  refer  the 
reader  to  in  [TR1],[TR2],[TR3],[TR4]  and  [TR5], 

In  addition  to  previous  reports,  in  this  work  we  also  analyze  the  problem  of  the  optimum 
sensor  selection  along  the  trajectory  of  a  tracked  target  in  a  multistatic  radar  systems.  To  this 
end  we  will  evaluate  the  Posterior  Cramer-Rao  Lower  Bound  (PCRLB)  of  the  sequential 
target  state  estimation  and  we  will  define  an  algorithm  that  exploits  this  mathematical  tool  to 
select  the  best  channels  for  the  tracking  of  the  target. 


-4- 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


2  The  Ambiguity  Function 


Classically,  the  properties  of  radar  waveforms  are  analyzed  and  presented  in  terms  of  the 
Complex  Ambiguity  Function,  originated  by  Woodward  in  the  1950s  [Wod80]. 

The  Complex  Ambiguity  Function  (CAF)  is  well  known  in  the  context  of  radar  as  a  key  tool 
for  determining  target  resolution  capability,  and  is  a  consequence  of  the  nature  of  the  optimal 
detector,  which  involves  decision-making  based  on  the  output  of  a  matched  filter  determined 
from  the  transmitted  waveform  [Tsa97].  As  a  matter  of  fact,  the  CAF  is  the  auto-correlation  of 
the  complex  envelope  of  the  waveform  with  a  copy  shifted  in  time  and  frequency,  and 
presents  the  point  target  response  of  the  waveform  as  a  two-dimensional  function  of  range  and 
Doppler,  showing  the  resolution,  sidelobe  structure  and  ambiguities  in  the  delay  and  Doppler 
domains.  The  CAF  is  intuitively  appealing  and  has  been  very  widely  used  -  indeed,  it  is  no 
exaggeration  to  say  that  every  serious  radar  engineer  on  the  planet  will  have  encountered  and 
used  the  Woodward  ambiguity  function. 

The  mathematical  definition  of  the  Complex  Ambiguity  Function  is  [Tsa97]: 

X  {XH  >*a’VH>Va)  =  jl  ~  Xa  V  (f  ~  XH  )  e~  ~K)' dt  (1) 


where  u(t)  is  the  complex  envelope  of  the  transmitted  signal,  xa  and  va  are  the  actual  delay  and 
Doppler  frequency  of  the  radar  target  respectively  and  xH  and  vH  are  the  hypothesized  delay 
and  frequency.  The  Ambiguity  Function  (AF)  is  defined  as  the  absolute  value  of  the  Complex 
Ambiguity  Function  and  is  clearly  maximum  for  xh  =  Ta  and  v//  =  va.  The  CAF  in  (1)  and  the 
AF  can  can  be  also  expressed  as  a  function  of  x  and  v,  where  x  =  xh  -  xa  and  v  =  vH  -  va.  In  this 
case  the  definition  of  the  Ambiguity  Function  is: 


A(t,v)  =  J  u(t)u (t-T)exp(-j27rvt)dt 


(2) 


Three  properties  of  the  AF  are  of  particular  interest  [Van71]. 

If  the  waveform  has  energy  E ,  then 

|-X’(t,v)|  <  |X(0,0)|  =  E  (3) 

Thus,  when  the  filter  is  matched  both  in  delay  and  Doppler  the  response  attains  the  maximum. 
If  the  filter  is  not  matched  then  the  response  assumes  a  value  lower  than  the  maximum. 

The  second  property  states  that  the  total  area  under  any  ambiguity  function  is  constant  and  it 
is  given  by 
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(4) 


J  J|X(t,v)|2  drdv  =  E2 

This  conservation  of  energy  statement  implies  that,  in  the  design  of  waveforms,  one  cannot 
remove  energy  from  one  portion  of  the  ambiguity  surface  without  placing  it  somewhere  else; 
it  can  only  be  moved  around  on  the  ambiguity  surface. 

The  third  property  is  a  symmetry  relation 

X(T,v)  =  X(-T-v)  (5) 

Moreover,  if  we  consider  the  CAF  for  v=0,  we  obtain  the  autocorrelation  function  of  u(t), 
similarly,  if  we  consider  the  CAF  for  r=0,  we  obtain  the  Fourier  transform  of  \u{t)\  . 

It  is  reasonable  to  ask  how  an  ideal  ambiguity  function  should  be.  The  answer  varies 
depending  on  the  aim  of  the  system  design,  but  a  common  goal  is  the  thumbtack  shape,  which 
features  a  single  central  peak,  with  the  remaining  energy  spread  uniformly  throughout  the 
delay-Doppler  plane.  The  lack  of  any  secondary  peak  implies  that  there  will  be  no  delay  or 
Doppler  ambiguities.  The  uniform  plateau  suggests  low  and  uniform  side  lobes,  minimizing 
target  masking  effects.  All  of  these  features  are  beneficial  for  a  system  designed  to  make  high 
resolution  measurement  of  targets  in  delay  and  Doppler,  or  to  perform  radar  imaging. 

In  our  work  we  evaluated  the  CAF  of  several  pulse.  In  particular  we  calculated  in  closed  form 
the  CAF  of  the  rectangular  pulse  [TR1],  the  burst  of  pulses  [TR1],  the  frequency/phase  coded 
waveform  [TR1],  the  Linear  Frequency  Modulated  (LFM)  pulse  [TR3]  and  the  Sinusoidal 
Frequency  Modulated  (SFM)  pulse.  For  more  detail  we  refer  the  reader  to  these  technical 
reports. 
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3  Frequency  Coded  Waveforms 

As  known,  in  order  to  have  an  auto-ambiguity  function  that  exhibits  a  narrow  thumbtack 
shape  with  low  sidelobes,  it  is  necessary  to  perform  a  pulse  compression  through  an  angle 
modulation.  In  this  Section  we  briefly  describe  the  Frequency-Coded  Waveforms  (FCW), 
where  the  transmitted  signal  is  frequency  modulated  through  a  code  sequence,  for  more 
details  we  refer  the  reader  to  [TR1].  The  complex  envelope  of  the  transmitted  signal  u(t)  is 
composed  of  N  pulses  p(t)  with  a  pulse  repetition  interval  Tr  and  it  can  be  written  in  the  form: 

N- 1 

u(t)  =  YJP(t-lTr)-e~j2*a,t  (6) 

/= o 

Where  f \=di)xc  and  {dN}={d0,dl,...,d,,...,dN_l}  is  a  sequence  of  integer  numbers  which 

belong  to  the  set  7a/={0,1,2,...,./V-1};  this  sequence  is  commonly  a  permutation  of  the  integer 
0,1,2,...  JV-l  and,  in  this  case,  the  code  sequence  is  a  full  code.  The  pulses  p(t)  are  often 
rectangular  pulses  with  time  duration  xc  and  pulse  repetition  interval  Tr  equal  to  xr.  In  this 
case,  the  transmitted  signal  u(t)  is  a  pulse  with  time  duration  xs=Nxc  composed  of  N  subpulses 
(or  chips)  with  time  duration  xc.  It  is  clear  that  an  FCW  signal  is  closely  related  to  the  code 
sequence  (or  placement  operators)  \d_N  } .  The  code  sequence  is  commonly  described  through 

the  geometric  array ,  that  is  a  matrix  whose  columns  represent  N  continuous  time  slices  and 
the  rows  represent  N  distinct  frequencies.  Through  the  geometric  array  it  is  possible  to 
calculate  the  sidelobes  matrix,  or  hits  matrix,  by  considering  the  geometric  array  as  a 
bidimensional  signal  and  calculating  its  autocorrelation  function  [TR1].  Clearly,  the  AF  of  an 
FCW  signal  is  closely  related  to  the  sidelobe  matrix  of  the  code,  in  fact  the  sidelobe  matrix 
describe  how  the  sidelobe  peaks  are  distributed  in  the  delay-Doppler  plane.  In  [TR1]  we 
described  the  techniques  commonly  used  to  construct  a  code  sequence  that  generates  a  good 
auto-ambiguity  function,  in  particular  we  analyzed  the  Costas-Sequences,  the  Pushing- 
Sequences,  the  Lee-Sequences,  the  Linear  Congruence  Codes,  the  Quadratic  Congruence 
Codes,  the  Cubic  Congruence  Codes  and  the  Hyperbolic  Congruence  Codes.  Anyway, 
particular  attention  has  been  devoted  to  the  analysis  of  the  techniques  to  construct  a  set  of 
code  sequence  that  generate  a  set  of  FCW  signals  for  netted  radar  system  application,  that  is  a 
set  of  signals  that  exhibit  a  good  auto-ambiguity  function  and  a  small  cross-correlation 
function  between  each  pair  of  signals  of  the  set.  In  particular  we  analyzed  the  Simulated- 
Annealing  technique  and  the  Extended  Quadratic  Congruences  (EQC)  technique,  that  will  be 
briefly  described  in  next  subSection.  For  more  detail  we  refer  the  reader  to  [TR1]. 
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3.1  Frequency  Coded  Waveforms  for  Netted  Radar  Systems 


With  the  ongoing  rapid  development  in  sensor  fusion  technology ,  netting  multiple  radar 
systems  through  fast  communication  links  is  becoming  a  trend  in  the  high-performance  radar 
design.  Via  information  fusion,  radar  performance  in  target  search,  tracking  and  recognition 
can  be  significantly  improved.  However,  the  waveforms  used  by  the  netted  radar  systems 
must  be  carefully  designed  to  avoid  the  self-interference  and  detection  confusion.  Then,  the 
waveforms  used  by  the  radars  in  a  complex  network  should  guarantee  an  auto-ambiguity 
function  that  exhibits  a  narrow  thumb  tack  shape  with  low  sidelobes  and,  in  contrast,  with 
small  cross-ambiguity  functions.  In  this  Section  we  describe  an  algebraic  technique  to 
construct  a  set  of  FCW  code  sequences  which  is  best  described  as  a  compromise  between  ease 

of  design  and  good  auto-  and  cross-ambiguity  properties.  The  set  of  sequences  {dn  are 
constructed  upon  an  extension  of  the  theory  of  quadratic  congruences  ( Extended  Quadratic 
Conguences  EQC)  [Bel88],  [Bel90],  [Bel91a],  [Bel91b],  [Tit81a],  [Tit81b],  [Tit91],  [Tit93], 
[Xia98].The  class  of  extended  quadratic  congruence  placement  operators  regroups  all  the 
operators  given  by  [TR1] 


dk  =  < 


k(k  + 1) 


0  <k< 


N- 1 


k(k  + 1)  .  ,  1 


-  +  (a -b)- 


N- 1 


<k<N- 1 


(7) 


2 

where  N  is  an  odd  prime  number  and  a  and  b  are  integers.  Note  also  that  the  identity  \N  ~lg=l 
holds  for  any  odd  integer,  so  that  (N  -l)/8  is  always  an  integer.  In  [TR1]  we  showed  that  the 

sequence  of  integers  {dn  is  a  permutation  of  the  set  { 1,2,. .  ,,N-l }  if  and  only  if  a  and  b  are 

not  both  quadratic  residues  ( QR )  or  quadratic  nonresidues  ( QNR )  of  the  odd  prime  N.  Each 
of  these  permutations  is  uniquely  defined  by  the  ordered  pair  (a,  b). 

The  number  X  is  a  QR  or  a  QNR  if 

U(A,-1)/2|  =(/ll  A0  =  {  1  QR  (8) 

1  'N  [-1  QNR 


where  (XIAO  is  the  so-called  Legendre  symbol. 

By  definition,  the  class  of  EQC  codes  regroups  all  code  words  parameterized  by  ordered 
pairs  (a,  b )  with  a  and  b  not  both  QR  or  QNR. 

As  there  are  (N- 1  )/2  QRs  and  (N- 1  )/2  QNRs,  there  are  exactly  (V- 1  y/2  EQC  code  words. 


-8- 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


As  showed  in  [TR1],  it  is  useful  to  organize  these  EQC  code  words  in  (N- 1  )/2  sets 
containing  N- 1  code  arrays  each,  with  the  property  that  each  unequal  pair  of  code  words 
among  the  N- 1  code  arrays  of  each  set  has  one,  and  only  one,  intersection  modulo  N  (at  k= 0 
by  construction).  In  general,  each  set  contains  the  ordered  pairs 
(aub]_),(bua2),(a2,b2),(b2,a3),---,(a(N-m,b(N-m)’(b(N-i)/2,ai)  for  some  sequences  of  all  distinct 
QR  or  QNR  a],  and  all  distinct  QNR  or  QR  b,  ;  in  addition,  when  (N-\)/2  is  odd,  one  set  can  be 
formed  with  the  ordered  pairs  (a\,N-a\),(N-a\,a\),. .  .,(a(N-\)i2,N-a(N-\)i2), 
( N-a(N.\)i2 ,  a(/V.|  ,/2).  As  an  example,  let  us  consider  the  case  N  =  5.  There  is  a  total  of  eight 
possible  EQC  code  words  that  can  be  generated  for  N  =  5,  which  are  grouped  in  two  sets 
containing  four  code  arrays  each.  The  first  set  is  shown  in  Figure  l.a;  it  corresponds  to  the 
sequence  of  pairs  (1,2),  (2,4),  (4,3),  and  (3,1),  respectively.  Note  that  there  is  only  one 
intersection  (at  k  =  0)  between  any  two  code  words  of  the  set.  Figure  1  .b  displays  the  second 
set  of  four  code  arrays;  it  is  associated  with  the  pairs  (1,3),  (3,4),  (4,2),  and  (2,1),  respectively. 

In  this  elementary  example,  the  only  alternative  sequence  is  obtained  by  simply  changing 
the  pair  ordering.  Again,  there  is  only  one  intersection  between  any  two  code  words  of  this 
set.  Any  set  different  from  those  of  Figure  1  would  not  possess  this  property. 

To  show  the  FCW-EQC  signals  guarantee  good  auto-  and  cross-ambiguity  properties,  Figure 
2-4  show  the  auto  and  cross  ambiguity  functions  of  the  FCW-EQC  signals  for  A=23  generated 

with  the  ordered  pair  (18,14)  and  (2,11),  namely  <^={0,  18,  8,  16,  19,  17,  10,  21,  4,  5,  1,  15, 
22,  20,  9,  12,  6,  14,  13,  3,  7,  2,  11}  e  dk={ 0,  2,  6,  12,  20,  7,  19,  10,  3,  21,  18,  17,  11,  16,  9, 
13,5,  8,22,  1,  14,  15,4}. 


Figure  1-  EQC  geometric  array  with  N=5. 
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Figure  2  -  Ambiguity  Function  of  sequence  dk . 
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Figure  3  -  Ambiguity  Function  of  sequence  dk . 
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Figure  4  -  Cross- Ambiguity  Function  of  sequences  dk  and  dk . 
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4  Adaptive  waveform  diversity  for  cross-channel  interference  mitigation 

In  previous  Section  we  showed  that  the  waveforms  used  by  the  radars  in  a  complex 
network  should  guarantee  good  target  detection  and  parameter  estimation  in  different 
scenarios  and  should  allow  an  optimal  access  to  the  same  transmit  channel.  We  also  showed 
that  FCW-EQC  waveforms  are  characterized  by  an  auto-ambiguity  function  that  exhibits  a 
narrow  thumb  tack  shape  with  low  sidelobes  and  a  good  cross-ambiguity  function,  almost  null 
in  the  delay-Doppler  plane. 

In  this  Section  we  analyze  a  scenario  composed  of  two  radars  transmitting  in  the  same 
frequency  band,  which  can  illuminate  the  same  area  looking  for  same  target.  These  two  radars 
can  use  either  the  same  or  different  codes.  In  particular,  we  investigate  the  impact  of  the 
presence  of  the  transmitted  signal  of  the  second  radar  (the  interfering  radar)  on  the  first  one 
(the  reference  radar)  in  the  estimation  of  the  target  Direction  Of  Arrival  (DOA).  The  analyzed 
scenario  is  pictorially  showed  in  Figure  5,  for  more  details  on  the  impact  of  the  intergfering 
signal  on  target  DOA  estimate  we  refer  the  reader  to  [TR1]  where  we  also  described  a  phased- 
array  FCW  receiver  and  we  evaluated  the  Cramer-Rao  Fower  Bounds  (CRFBs)  of  target 
DOA  estimate. 


Figure  5  -  Impact  of  the  interfering  radar  on  the  reference  radar. 
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In  a  typical  phased  array  radar,  a  single  beam  is  formed  on  transmission  and  two  or  more 
beams  are  formed  on  reception.  We  assume  here  that  the  system  is  a  linear  array  radar  that 
estimates  the  target  DOA  by  using  the  sum  channel  2  on  transmission  and  two  matched 
channels,  the  sum  2  and  the  difference  A  on  reception.  The  two  channels  fz(dT)  and  f\(0T),  or 
antenna  patterns,  are  defined  as  the  complex  amplitude  profiles  versus  target  azimuth  angle  9t 
(for  more  detail  we  refer  the  reader  to  [TR1]  and  [TR2]).  We  assume  that  the  sum  channel  has 
a  HPBW  of  3°  and  that  the  reference  system  is  a  tracking  system,  then  the  target  and  the  radar 
are  static  during  the  Time  on  Target  (ToT).  If  the  radar  use  a  Pseudo-Monopulse  (PM) 
estimator  as  in  [TR1],  it  transmits  a  sequence  of  K  pulses  used  for  DOA  estimation. 

The  performance  of  PM-DOA  estimator  in  presence  of  an  interfering  radar  has  been  analyzed 
in  detail  in  [TR1],  where  we  also  analyzed  the  performance  of  the  ML  (Maximum  Likelihood) 
estimator,  not  reported  here  for  the  lack  of  space. 

As  showed  in  [TR1],  the  interfering  radar  introduces  on  the  DOA  estimation  a  strong  bias  that 
depends  on  the  angular  position  of  the  interference  and,  particularly,  on  the  ratio  I  C/A  I,  where 
C  is  the  complex  value  of  Cross-Ambiguity  Function  at  the  sampling  instant  while  A  is  the 
complex  value  of  Auto- Ambiguity  Function  at  the  sampling  instant. 

In  [TR2]  we  proposed  a  technique  for  mitigating  the  effect  of  the  interfering  radar  that  will  be 
briefly  described  in  this  Section.  For  more  details  we  refer  the  reader  to  [TR2]. 

The  proposed  idea  is  simple,  that  is  if  the  target  is  still  during  the  ToT,  and  there  is  not 
interference,  the  difference  between  two  subsequent  DOA  estimates  is  related  only  to  the 
power  of  disturbance  (noise  and  clutter)  contribution  and  for  high  signal-to-disturbance 
power  ratio  tends  to  zero.  If  there  is  an  interference  and  the  contribution  of  the  interfering 
signal  has  a  different  weight  on  each  estimate  the  difference  between  the  two  subsequent 
DOA  estimates  depends  also  on  the  interference. 

Then  the  difference  between  two  estimates  can  be  used  as  key  statistic  to  detect  the  presence 
of  the  interference  and  waveform  diversity  in  transmission  should  be  applied  for  having 
different  contributions  of  the  interfering  signal  on  each  estimate.  For  this  purpose,  during  the 
ToT  the  radar  transmits  two  bursts  of  N  frequency  hop  coded  pulses.  The  frequency  hop 
signals  used  in  the  two  bursts  are  generated  with  different  codes.  The  receiver  is  initially 
matched  to  the  first  frequency  hop  pulse  used  in  the  first  burst  and,  in  a  second  time,  it  is 
matched  to  the  second  frequency  hop  pulse  used  in  the  second  burst.  Based  on  this,  it  is 
possible  to  write  the  expression  of  the  two  pairs  of  signals  received  on  the  sum  and  on  the 
difference  channels.  Let’s  suppose  that  the  Doppler  frequencies  of  the  target  and  of  the 
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interference  radar  are  the  same  with  respect  to  the  reference  radar.  In  this  case,  the  reference 
radar  cannot  distinguish  between  the  two  signals  by  means  of  Doppler  processing. 

Without  lack  of  generality,  we  assume  that  both  Doppler  frequencies  are  zero;  then  the 


received  signals  can  be  written  as: 

We (n)  =  a-  A ■  (4 )  +  p ■  C,  ■  fL (&, )  +  dxx (n) ,  (9) 

y1A(n)  =  a-A-fz(tfT)-fA(&T)  +  /3-Cl-fA(tfI)+dlA(n),  (10) 

y2x(n)  =  a' A2(4)  +  Pci-  fM)  +  d2X(n)  (li) 

W.A (»)  =  a ■  A  ■  /e (4 )  •  /A  (4 )  +  p  •  c2  ■  fA  )  +  d2  A ( n )  (12) 


Subscripts  “1”  and  “2”  refer  to  the  first  and  the  second  burst,  which  are  received  in  series, 
and  n=0,l,...)V-l.  The  first  term  in  y,-,x  and  v,.a  (/=  1 ,2)  is  due  to  the  target  signal;  it  depends  on 
the  target  DOA  0T  through  f  y(0T)  in  sum  channel  and  through  fy{0r)f\(0r)  in  the  difference 
channel.  This  is  due  to  the  two-way  antenna  gain.  The  second  term  is  due  to  the  signal 
transmitted  by  the  second  radar,  which  has  DOA  0j ,  and  it  depends  on  the  one-way  gain  of  the 
antenna  patterns,  i.e.  fz(di)  and  fA/)/)  fA ( f), ) .  dix  and  c/,.A  are  the  noises  on  the  two  channels, 
modeled  as  mutually  independent  complex  white  Gaussian  processes.  A  is  the  complex  value 
of  the  Auto- Ambiguity  Function  of  the  transmitted  signal,  C i  is  the  complex  value  of  the 
Cross-Ambiguity  Function  between  the  first  frequency  hop  pulse  and  the  interfering  signal 
and  Ct  is  the  complex  value  of  the  Cross-Ambiguity  Function  between  the  second  frequency 
hop  pulse  and  the  interfering  signal.  If  the  reference  radar  is  synchronized  in  reception  and 
transmission  and  the  receiver  is  tuned  on  the  Doppler  of  the  target,  then  A  is  the  energy  of  the 
transmitted  pulse;  a  is  the  complex  amplitude  of  the  received  echo  and  /?  is  the  complex 
amplitude  of  the  interfering  signal.  They  are  modeled  as  complex  Gaussian  mutually 
independent  random  variables;  in  short  a~CM0,cr(<J  and  fi~CN(0,o~p).  The  signal-to-noise 
power  ratio  is  defined  as  SNR=\A\2o2a/o2d  and  the  signal-to-interference  ratio  as  SIR=o2Ja2i g. 

Then,  using  the  N  samples  received  on  each  of  the  two  channels  2  and  A  from  the  first 
burst,  the  signal  processor  performs  a  first  target  DOA  estimation  through  the  Pseudo- 
Monopulse  (PM)  technique.  In  particular,  the  signal  processor  forms  the  monopulse  ratio 
defined  by  r\(n)=Re {y / tA(n)/y / x(n) }  for  each  pulse,  where  Re{}  denotes  the  real  part.  In 
absence  of  disturbance  and  in  the  presence  of  only  one  target,  the  monopulse  ratio  reduces  to 
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ri(n)=Re{fA{OT)/fz.(OT)}  from  which  the  angular  location  of  the  target  for  each  pulse  can  be 
determined.  Finally,  we  obtain  ^  =  r:  1  / N ,  where  is  the  target  DOA  estimate  for 

each  pulse.  Similarly,  the  signal  processor  performs  a  second  target  DOA  estimation  z?2  using 

the  N  samples  in  Eqs.  (11)-(12)  received  on  the  two  channels  from  the  second  burst.  The 
proposed  Modified  Pseudo-Monopulse  (MPM)  estimate  is  obtained  from  these  two  PM 
estimates  as  described  in  the  following  (see  [TR2]  for  more  details).  First  the  difference 

AiX  =  -  $2  is  calculated  and  then  compared  with  a  threshold  X.  If  A0<A  the  MPM  estimate 

/V  /V  A  /  A  A  \  / 

is  calculated  as  the  average  of  the  two  PM  estimates,  ^  and  ,  that  is  $7  .  =  I  ^  + 1\  )  2 . 

On  the  contrary,  if  A0>a,  a  third  PM  estimate  is  calculated  by  using  the  samples 
received  on  the  two  channel  E  and  A  from  a  third  burst  of  N  frequency  hop  pulses,  coded  whit 
a  third  code  sequence.  Then,  the  signal  processor  calculates  A&{=  and 

A  x\  =  .  If  the  lowest  between  these  two  quantities  is  less  than  X,  then  the  MPM 

estimation  of  the  target  DOA  is  calculated  as 

A  A 

t9  +  7 9 

rtVfi  Atf<AA 

A=  .  2.  "  (13) 

A^>A^2 

^  2 

If  both  A0i  and  A02  are  greater  than  X,  the  MPM  estimation  of  the  target  DOA  is  calculated 

A  /  A  A  A  \  / 

as  the  average  of  the  three  PM  estimates,  that  is  i\  =  ( + 1%  +  ?2, 1  /3  . 

In  this  last  case,  the  target  DOA  can  be  also  calculated  as  the  value  closest  to  among 
the  three  PM  estimates,  where  &p  is  the  MPM  estimation  of  the  target  DOA  performed  in  the 
previous  tracking  time  lag.  However,  it  is  important  to  consider  that  the  value  of  $p  is  not 

always  available  and  it  is  often  difficult  to  perform  the  initialization  of  . 

To  evaluate  the  impact  of  the  presence  of  the  interfering  radar,  the  RMSE  of  the  DOA 
estimator  has  been  derived  by  running  104  Monte  Carlo  simulations.  In  our  simulations,  we 
set  ()/—0°,  N=  8,  SIR=0dB  and  SNR= 20dB.  With  this  value  of  SNR  the  probability  of  detection 
of  the  filter  matched  to  the  transmitted  signal  is  almost  unitary,  even  for  Pfa=  10'6.  Due  to  the 
high  value  of  SNR,  the  performance  of  the  DOA  estimator  are  mainly  affected  by  the  presence 
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of  the  interfering  radar.  Let’s  first  show,  as  comparison  for  the  MPM  estimator,  some  results 
relating  to  the  PM  estimator  in  presence  of  interference.  In  the  case  of  PM  estimation  the 
received  signals  are  [TR2] : 

yL(n)  =  aAfz(&T)  +  /3C  fM)  +  dz(n)  (14) 

yA{n)  =  a  A  /e(4)./a(4  )  +  Pc  fM )  +  dA(n)  (15) 

In  our  analysis,  we  considered  0<IC/AI<1.  The  worst  case  is  when  IC/4I=1;  this  value 
characterizes  two  synchronized  radars  using  the  same  code.  The  best  case  is  when  IC/AI=0; 
this  value  characterizes  the  case  of  synchronized  radars  using  two  orthogonal  codes.  In  Figure 
6  we  plotted  the  RMSE  of  the  PM  estimator  as  a  function  of  the  interfering  DOA.  We  suppose 
that  the  PM  estimator  performs  the  estimation  of  the  target  DOA  using  a  burst  of  24 
frequency  hop  pulses  (that  corresponds  to  three  bursts  of  N=  8  pulses).  The  RMSE  has  been 
calculated  for  IC/AI=1  and  IC/AI~I/(0,1),  that  is  when  IC/AI  is  an  uniformly  distributed  random 
variable  in  the  interval  [0,1].  As  apparent  in  Figure  6,  the  performances  of  the  PM  estimator 
are  strongly  dependent  on  the  interfering  signal. 


0  1  2  3  4  5  6 

8  [degrees] 

Figure  6  -  RMSE  of  the  PM  estimator  as  a  function  of  the  interfering  DOA. 
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In  Figure  7  we  reported  the  RMSE  of  the  MPM  estimator  as  a  function  of  the  interfering 
DOA  for  different  values  of  the  threshold  1.  In  this  simulation,  when  both  A0i  and  A02  are 
greater  than  X,  the  MPM  estimation  of  the  target  DOA  has  been  calculated  as  the  value  closest 

to  z?  among  the  three  PM  estimates.  The  value  of  z?  has  been  modeled  as  a  Gaussian 

random  variables  with  mean  equal  to  the  target  DOA  and  variance  equal  to  the  Cramer-Rao 
Lower  Bound  (CRLB)  of  the  target  DOA  estimation  performed  using  a  burst  of  2N  pulses  (for 
more  details  see  [TR2]).  In  this  simulation  we  also  set  ICi/AI=l,  IC2/AI=0  and  I  CyA  1=0.  That 
is,  we  supposed  that  the  codes  transmitted  in  the  first  burst  by  the  reference  radar  and  the 
interfering  radar  are  the  same,  but  the  codes  transmitted  in  the  second  and  in  the  third  burst 
are  orthogonal  to  the  interfering  signal.  As  apparent  in  Figure  7,  when  the  threshold  X  is  large 
(i.e.  X>2)  the  performance  of  the  MPM  estimator  is  almost  the  same  as  the  PM  estimator.  This 
is  due  to  the  fact  that  in  this  case  the  probability  that  A9  is  lower  than  X  Pr{  A0<1}  is  almost 
unitary  and  the  MPM  estimator  works  as  a  PM  estimator  which  performs  the  target  DOA 
estimation  using  a  burst  of  2 N  pulses.  The  only  difference  is  that  the  MPM  estimator  uses  two 
bursts  of  N  frequency  hop  pulses  coded  with  two  different  code  sequences  instead  of  a  burst 
of  2 N  frequency  hop  pulses  coded  with  the  same  sequence. 

On  the  other  hand,  when  the  threshold  X  decreases,  the  probability  Pr{A0<l}  decreases  and 
the  MPM  estimator  calculates  the  target  DOA  using  two  or  three  bursts  of  N  pulses.  It  is 
apparent  from  Figure  7  that  the  lower  is  X,  the  lower  is  the  dependence  of  the  estimation 
accuracy  on  the  interfering  DOA.  When  the  threshold  X  tends  to  zero,  also  the  probability 
Pr{  A0<1}  tends  to  zero,  and  the  MPM  estimator  always  processes  three  different  PM  DOA 
estimates.  It  is  clear  from  Figure  7  that  this  algorithm  minimizes  the  RMSE,  that  is,  when  the 

estimate  d-p  is  available,  the  optimum  value  of  X  is  0.  However,  the  value  of  d-p  is  not  always 

available,  then,  in  the  more  general  case,  when  both  A0i  and  A92  are  greater  than  X,  the  MPM 
estimation  of  the  target  DOA  must  be  calculated  using  the  three  PM  estimates. 

In  Figure  8  we  report  the  RMSE  of  the  MPM  estimator  as  a  function  of  the  interfering  DOA 
for  different  values  of  the  threshold  X.  In  this  simulation,  when  both  A0i  and  A92  are  greater 
than  X,  the  MPM  estimation  of  the  target  DOA  has  been  calculated  as  the  average  of  the  three 

PM  estimates,  that  is  z7r  =  (z?,  +  z?2  +  ?73  .  Also  in  this  simulation  we  set  \C\/A\=\,  !C2/AI=0 

and  IC3/AI=0.  In  this  case,  when  1=0,  the  performance  of  the  MPM  estimator  is  strongly 
dependent  on  the  interfering  signal.  Through  simulations,  it  is  possible  to  find  the  value  of  X 
that  minimizes  the  RMSE,  in  our  case,  as  apparent  from  Figure  8,  this  value  is  1=0.3. 
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Then,  using  the  MPM  estimator  with  a  proper  threshold,  it  is  possible  to  improve  greatly  the 
performance  of  the  radar  with  respect  to  the  PM  DOA  estimator. 

The  case  of  iC2/AI=0  and  IC3/AI=0  is  the  best  for  the  second  and  third  bursts.  In  Figure  9  we 
show  some  results  with  more  general  cases.  In  this  Figure  the  RMSE  of  the  MPM  estimator 
has  been  calculated  by  assuming  ICi/AI=l,  IC2/AI~t/(0,l-y)  and  IC3/AI~t/(0,l-y).  The  results 
are  reported  as  a  function  of  the  interfering  DOA,  for  different  values  of  the  parameter  y  and 
for  ?i=0.3.  Also  in  this  simulation,  when  both  A0i  and  A02  are  greater  than  X,  the  MPM 
estimation  of  the  target  DOA  has  been  calculated  as  the  average  of  the  three  PM  estimates. 

It  is  apparent  that,  for  high  values  of  y  the  MPM  estimator  allows  again  an  improvement  with 
respect  to  the  MP  estimator.  It  is  important  to  notice  that,  using  proper  frequency  hop  codes, 
the  values  of  I  C/A  I  tends  to  zero  for  every  Doppler-delay  shift  between  the  reference  and  the 
interfering  signals.  Then  the  value  of  y  commonly  tends  to  one  and  the  results  obtained  for 
low  values  of  y  have  to  be  considered  the  worst  cases. 
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Figure  7  -  RMSE  of  the  MPM  estimator  as  function  of  0/  for  different  value  of  L 


Figure  8  -  RMSE  of  the  MPM  estimator  as  function  of  0/  for  different  value  of  L 


Figure  9  -  RMSE  of  the  MPM  estimator  as  function  of  0/  for  different  value  of  y. 
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5  Bistatic  Radar  Systems  and  Bistatic  Ambiguity  Function 

In  previous  Section  we  analyzed  a  netted  system  composed  by  monostatic  radar  system. 
In  the  most  general  case,  in  a  multistatic  radar  system  the  receivers  can  also  exploit  the  signal 
emitted  by  non  colocated  transmitters.  To  this  end,  it  is  very  important  to  define  and  analyze 
bistatic  radar  system.  This  was  the  main  topic  of  [TR3]  and  [TR4]  summarized  hereafter. 
Bistatic  radar  may  be  defined  as  a  radar  in  which  the  transmitter  and  receiver  are  at  separate 
locations.  The  very  first  radars  were  bistatic,  until  pulsed  waveforms  and  T/R  switches  were 
developed.  Since  then  interest  has  varied  up  and  down,  but  is  demonstrably  now  at  a  high 
level,  with  numerous  experimental  systems  being  built  and  the  results  reported  in  the 
literature.  Rather  fewer  operational  systems,  though,  have  been  deployed. 

Bistatic  radars  can  operate  with  their  own  dedicated  transmitters,  which  are  specially  designed 
for  bistatic  operation,  or  with  transmitters  of  opportunity,  which  are  designed  for  other 
purposes  but  found  suitable  for  bistatic  operation.  When  the  transmitter  of  opportunity  is  from 
a  monostatic  radar  the  bistatic  radar  is  often  called  a  hitchhiker.  When  the  transmitter  of 
opportunity  is  from  a  non-radar  transmission,  such  as  broadcast,  communications  or 
radionavigation  signal,  the  bistatic  radar  is  called  Passive  Coherent  Location  (PCL). 

Finally,  transmitters  of  opportunity  in  military  scenarios  can  be  designated  either  cooperative 
or  non-cooperative,  where  cooperative  denotes  an  allied  or  friendly  transmitter  and  non- 
cooperative  denotes  a  hostile  or  neutral  transmitter.  Passive  bistatic  radar  operations  are  more 
restricted  when  using  the  latter.  In  this  Section  we  first  define  the  bistatic  parameters  and  the 
bistatic  coordinate  system  and  then  we  analyze  the  Bistatic  Ambiguity  Function. 


5.1  Bistatic  Geometry 

Before  starting,  it  is  necessary  to  describe  the  coordinate  system  used  to  represent  a 
bistatic  radar  geometry.  Figure  10  shows  the  coordinate  system  and  its  parameters.  The 
positions  of  the  TX,  RX  and  of  the  target  are  generic.  Considering  an  ordinary  Cartesian  grid, 
the  TX  is  located  at  point  T,  whose  coordinates  are  (xT,  y-/j,  the  RX  is  located  at  point  R  in 
(. xr,  yR)  and  the  target  is  located  at  point  B,  whose  coordinates  are  (x,  y).  The  triangle  formed 
by  the  transmitter,  the  receiver  and  the  target  is  called  the  bistatic  triangle. 

As  shown  in  Figure  10,  the  sides  of  the  bistatic  triangle  are  RT,  Rr  and  L,  where  R/  is  the 
range  from  transmitter  to  target,  Rr  is  the  range  from  receiver  to  target  and  L  is  the  baseline 
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between  the  transmitter  and  the  receiver.  The  internal  angles  of  the  bistatic  triangle,  that, 
without  lack  of  generality,  are  assumed  to  be  positive,  are  a,  f>  and  y.  In  particular,  the  bistatic 
angle  [>  is  the  angle  at  the  apex  of  the  bistatic  triangle,  at  the  vertex  which  represents  the 
target.  Assuming  that  the  coordinates  of  the  transmitter,  the  receiver  and  the  target  are  known, 
it  is  possible  to  calculate  all  the  parameters  of  the  bistatic  triangle.  0/  and  Or  are  the  look  angle 
of  the  transmitter  and  the  look  angle  of  the  receiver,  respectively,  they  are  measured  positive 
clockwise  from  the  vector  normal  to  the  baseline  pointing  towards  the  target. 

From  Figure  9.1,  we  have  that  dT=  90°-a,  Or  ■  y-90°,  ft  =  l80°-a-y=0j—0R,  and  from  the 
cosine  law  we  obtain:  Rj  =R  if +L~ +2R  aLfmO  r.  which  gives  the  range  from  transmitter  to 
target  Rt,  as  a  function  of  the  range  from  receiver  to  target  Rr  and  the  look  angle  of  the 

receiver  Or.  Figure  9.1  also  shows  the  target  velocity  vector  V  ;  (j)  is  the  angle  between  the 
target  velocity  vector  and  the  bistatic  bisector,  which  is  measured  in  a  positive  clockwise 
direction  from  the  bisector.  In  particular  the  bistatic  bisector  is  represented  by  the  vector  Bl , 
where  1  is  the  incenter  of  the  bistatic  triangle,  whose  coordinates  are  ( xj ,  yj). 

The  coordinates  of  the  incenter  can  be  easily  obtained  as 


(■ */>?/ ): 


L  +  Rr  +  Rj 


(*.?)+ 


L  +  Rr  +  Rt 


(^7 ’’  Jr  ) 


L  +  Rr  +  Rj 


(-W*) 


Figure  10  -  Bistatic  geometry. 
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In  the  bistatic  geometry,  an  important  parameter  is  the  radial  velocity  Va,  which  is  the  target 
velocity  component  along  the  bistatic  bisector.  From  the  observation  of  Figure  9.1,  we  obtain 

Vicos </> .  Using  the  notation  V  =  Vx  ■  x  +  Vy  ■  y  ,  it  is  easy  to  verify  that: 


V  =VBI/ 


Bl 


(■v/-.r)yt  +  (y/-y)y>, 

yl(x,-x)2  +  (y,-y)2 


(17) 


The  bistatic  radar  geometry  can  be  completely  specified  in  terms  of  any  three  of  the  five 
parameters,  0T,  0R,  L,  Rr  and  R T.  In  this  Section  we  will  use  0R,  L,  and  Rr  that  can  be  obtained 
using  the  following  equations: 


Z,  \jixj-  xR )  +  ( yT  y  R ) 

rr  =^(-x~Xr)2  +(y-yR)2 

Rt  =yl(x-xT)2  +(y-yT)2 


&R  =  cos" 


f  Rl+L2-R2^ 

V  2RrL  J 


71 

~2 


(18) 

(19) 

(20) 

(21) 


5.2  Bistatic  Ambiguity  Function 

As  showed  in  Section  1,  the  mathematical  definition  of  the  Ambiguity  Function  is: 


X(rH,ra,VH, va  )|  =  \+lu(t -  ra )u  (t-TH) cxp(-j27r(vH  -va)t)dt 


(22) 


where  u{t)  is  the  complex  envelope  of  the  transmitted  signal,  za  and  va  are  the  actual  delay  and 
Doppler  frequency  of  the  radar  target  respectively  and  xH  and  vH  are  the  hypothesized  delay 
and  frequency.  The  AF  in  (22)  can  be  also  expressed  as  a  function  of  x  and  v,  where 
x  =  xH  -  Ta  and  v  =  vH  -  va.  As  known,  in  the  monostatic  case  there  is  a  linear  relationship 
between  xa  and  va,  and  the  range  position  Ra  and  radial  velocity  Va  of  the  target,  more 
specifically  xa  =2 RJc  and  v„  =-2Vclfc/c.  Similar  relations  hold  for  xh  and  vH.  Due  to  this  linear 
relationship,  the  AF  in  the  range  -  velocity  plane  has  the  same  behavior  of  the  one  expressed 
as  a  function  of  x  and  v,  except  for  a  scale  factor.  Therefore,  in  the  monostatic  configuration, 
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the  information  about  the  target  delay  and  the  target  Doppler  shift  directly  provides 
information  about  the  target  range  and  the  target  velocity.  This  is  different  in  the  bistatic  case, 
where  the  relation  between  time  delay  and  Doppler  frequency,  and  target  distance  and 
velocity  is  not  linear.  Referring  to  the  bistatic  geometry  of  Figure  10,  for  obtain  the 
expression  of  the  bistatic  ambiguity  function,  we  must  replace  in  (22)  the  relations  [Tsa97]: 


th(Rr,Or,L) 


Rr  +  \]  Rp>  +  L  +2RRLsinOR 

c 


(23) 


So 


1 


vh(Rr,Vb,0r,L)  =  2^VbI-  + 


Rr  +  L  sin  0R 


2  2^  Rr+  1}  +  2RSL  sin  0R 


(24) 


Similar  relations  hold  for  th  and  vh  ■  It  is  clearly  apparent  from  equations  (23)  and  (24)  that  in 
the  bistatic  case,  the  Doppler  shift  and  the  delay  depends  on  the  geometry  of  the  bistatic 
triangle  and  the  relation  between  time  delay  and  Doppler  frequency,  and  target  distance  and 
velocity  is  not  linear.  Due  to  these  non  linear  equations,  it  is  apparent  that  the  Bistatic 
Ambiguity  Function  depends  also  on  the  bistatic  geometry  parameters,  i.e.  the  target  direction 
of  arrival,  the  bistatic  baseline  length  and  the  distance  between  the  target  and  the  receiver. 
This  dependence  is  very  strong  and  can  be  appreciated  later  with  an  illustrative  example, 
considering  the  case  of  a  target  close  to  the  baseline  joining  the  transmitter  and  receiver. 

In  [TR3]  we  analyzed  the  performance  when  the  transmitted  waveform  is  a  burst  of  Linear 
Frequency  Modulated  (LFM)  pulses,  while  in  [TR4]  we  also  analyzed  the  performance  when 
the  transmitted  waveform  is  a  Sinusoidal  Frequency  Modulated  (SFM)  pulse.  Hereafter  we 
will  consider  the  case  of  a  burst  of  LFM  pulses,  for  the  SFM  case  we  refer  the  reader  to 
[TR4].  As  calculated  in  [TR3],  the  Complex  Ambiguity  Function  (CAF)  of  a  burst  of  LFM 
pulses  can  be  expressed  as: 

1  ,  sinf^v(iV-|/x|)r„) 

X{t.v)  =  -  £  exp (jzv(N-l+n)Ts)  \)  'x,(T-„r,v)  (25) 

N  nJw- 1)  sin  \7ZVTr  ) 


where 


Xl{r,v)  = 


sin  (^T(v-A:r)(l-|r|/r)) 
7dT{v-kr) 


red 


f  T  3 

2  T 


(26) 


is  the  Complex  Ambiguity  Function  of  a  single  LFM  pulse. 


1  Note  that  Ra  and  Va  are  the  actual  range  and  bistatic  velocity,  while  RT  and  Vl;=Vcos(p  are  the  hypothesized 
range  and  bistatic  velocity 
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Moreover,  if  we  limit  the  delay  to  the  mainlobe  area,  namely  to  Irl  <  T,  the  AF  reduces  to: 


X(r.v) 


sin(^r(v-kr)(l-  t/t)) 

sin  {tzvNTr) 

7fr(v-kr) 

N  sm(7rvTR) 

for  \t\  <  T 


(27) 


To  link  this  equation  to  the  bistatic  geometry  of  Figure  10  and  to  obtain  the  expression  of  the 
Bistatic  AF  (BAF),  we  must  replace  in  (27)  the  relations  z  =  zh  -  za  and  v  =  vh  -  va  calculated 
using  (23)  and  (24).  The  contour  plot  of  the  BAF  is  illustrated  in  Figure  11  in  the  plane  Rr-Vb 
with  VB=Vcoscp,  14=600 m/s  and  Ra= 20Km  and  L=50Km.  The  presence  of  discrete  peaks 
(nails)  is  evident  even  in  the  bistatic  plane,  even  if  they  are  not  symmetrically  distributed. 

The  main  peak  corresponds  to  14=600m/s  and  /?a=20Km.  The  shape  of  the  bistatic  function 
strongly  depends  on  the  target  angle  Or.  particularly  for  high  values  of  BT.  To  highlight  this 
phenomenon,  in  Figures  12-15  we  show  the  zero-delay  and  zero-Doppler  slices  of  the 
ambiguity  function  for  different  values  of  9r  and  BT.  Both  cuts  are  maximum  for  values  of 
range  and  target  velocity  corresponding  to  the  true  values  14=600m/s  and  /4=20Km. 

For  values  of  0R  close  to  — zr/2  the  bistatic  AF  presents  multiple  peaks.  The  worst  case  is  for 
0r=-tz/2,  that  is,  when  the  target  is  on  the  baseline.  If  the  target  is  between  the  transmitter  and 
the  receiver,  the  AF  is  flat  and  the  range  and  velocity  resolutions  are  completely  lost.  For 
values  of  Or  far  from  -nil  the  shape  of  the  bistatic  ambiguity  function  is  practically  the  same 
(see,  for  instance,  Or  =  -tc,  and  Or  =  n/6  in  Figs.  12-15).  For  increasing  values  of  N  the  range 
resolution  improves,  but  many  peaks  appear  in  the  bistatic  AF  shape. 
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Figure  11  -  BAF  of  a  burst  of  Linear  Frequency  Modulated  pulses,  7? 7=250,  T^lms, 
7=250ps,  N= 8,  47 n  ,  L=50Km,  Ffl=600m/s,  7?fl=20Km. 


Figure  12  -  BAF  of  a  burst  of  Linear  Frequency  Modulated  pulses;  Zero-Doppler  cut, 
BT= 20,  7>lms,  7=250ps,  N=  8,  L=50Km,  Fa=600m/s,  7?fl=20Km. 
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Figure  13  -  BAF  of  a  burst  of  Linear  Frequency  Modulated  pulses;  Zero-delay  cut, 
BT= 20,  Tr-  1  ms,  T=250ps,  N= 8,  L=50Km,  Ffl=600m/s,  fla=20Km. 


Figure  14  -  BAF  of  a  burst  of  Linear  Frequency  Modulated  pulses;  Zero-Doppler  cut, 
BT= 250,  7  /^  lms,  T=250ps,  N= 8,  L=50Km,  Ffl=600m/s,  7?fl=20Km. 
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Figure  15  -  BAF  of  a  burst  of  Linear  Frequency  Modulated  pulses;  Zero-delay  cut, 
BT= 250,  7  /^lms,  F=250ps,  N=  8,  L=50Km,  Fa=600m/s,  Ra= 20Km. 
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6  Channel  Performance  Evaluation  in  a  Bistatic  Radar  System 


6.1  Bistatic  Cramer-Rao  Lower  Bounds 


The  AF  directly  determines  the  capability  of  a  system  to  resolve  two  targets  that  exist  at 
different  ranges  from  the  radar  and  have  different  radial  velocities.  When  the  receiver  signals 
from  the  target  have  similar  energy,  the  resolution  is  equal  to  the  half  power  width  of  the  AF 
mainlobe.  The  AF  is  also  related  to  the  accuracy  with  which  the  range  and  the  velocity  of  a 
given  target  can  be  estimated.  When  the  Signal  to  Noise  power  Ratio  (SNR)  is  high,  the 
CRLBs  on  estimation  accuracy  are  dependent  on  both  the  SNR  and  the  second  derivatives  of 
the  AF,  that  is,  the  sharpness  of  the  AF  mainlobe.  Unlike  the  ambiguity  function  which 
provides  information  on  the  global  resolution,  the  CRLBs  are  a  local  measure  of  estimation 
accuracy.  Anyway,  both  can  be  used  to  asses  the  error  properties  of  the  estimates  of  the  signal 
parameters.  In  [Van71]  the  author  derived  a  relationship  between  CRLB  and  ambiguity 
function,  which  has  been  successfully  used  in  the  analysis  of  passive  and  active  arrays 
[DogOl].  In  the  monostatic  configuration,  [Van71]  claims  that  for  the  Fisher  Information 
Matrix  (FIM)  the  following  relationship  holds  (for  more  details  see  [TR3]): 


Jm{t,v)  =  -2SNR 


d20(r,v) 

d20(r,v) 

dr2 

dtdv 

d2@(r,v) 

d20(r,v) 

dvdr 

dv2  J 

=  -2SNR  J 


(28) 


T= 0,K=0 


I  |2 

where  ®(t,v)  =  \X(t,v)\  and  SNR  is  the  signal-to-noise  power  ratio  at  the  receiver.  The  AF 

is  the  heart  of  this  expression  since  it  is  the  log-likelihood  function  excluding  the  effect  of 
signal  attenuation  and  clutter.  In  [TR3]  we  report  the  proof  of  relation  (28).  The  property  in 
(28)  does  not  depend  on  the  choice  of  the  parameters  of  the  ambiguity  function,  then  it  holds 
for  both  monostatic  and  bistatic  case.  From  (28)  the  CRLBs  follow: 

CRLB(rJ  =  [J„fe,r,,)]-;  and  CRLB(n„)  =  [J  m (TaNa)\7\-  In  the  bistatic  configuration  we 

should  write  the  ambiguity  function  in  terms  of  the  bistatic  t(Rr,9r,L )  and  v(Rr,  Vb,Or,L )  and 
derive  it  with  respect  to  the  useful  parameters  Rr  and  VR.  Then 
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J,(Rg.VB)  =  -2SNR 


(29) 


>e(R„v„) 

z'MRM] 

3R/ 

y@(RK.vt) 

32®(«*V.) 

5Vb3Rk 

j 

For  the  calculation  of  the  CRLBs  in  the  bistatic  domain  we  can  partially  use  the  results  of  the 
monostatic  domain.  Following  the  “chain  rule”  (see  [TR3]  for  details)  of  the  derivative  we 
can  prove  that: 

j,(/;s,v',)=pj„(r,K)pr  (30) 


where 


dr 

dv 

3RS 

3Rr 

dr 

dv 

and 


_  fc 

dRp  2c 


V„ 


L  cos  0„ 


^  Rr  +  L~  +  2RxLsin  0R  j 


3/2  1 


-  + 


Rr  +  L  sin  0 R 


3  2yj  Rr  +  L  +  2  RrL  sin  0R 


(31) 


(32) 


dv  =  2./;.  1 1  |  Z^+Lsin^ 

c  p  2^R2  +  L2  +  2RRLsin0R 


dr  1 

f 

1  +  _ 

Rr+L  sin  0R 

dRR  c 

±  l 

v  1 

J Rr  +  L  +  2 RrL  sin  0R  ^ 

dr 

dVB 


=  0 


(33) 


(34) 


(35) 


The  relationship  showed  in  (30)  is  very  interesting  because  the  two  effects  that  describe  the 
Bistatic  FIM  are  separated.  In  particular,  matrix  P  takes  into  account  only  the  effect  of  the 
bistatic  geometry  while  J m(fv)  takes  into  account  only  the  effect  of  the  transmitted 
waveform.  The  Cramer-Rao  lower  bounds  are  given  by  the  inverse  of  the  Fisher  Information 
Matrix,  therefore 


CRLB(RS)  = 


[J»L 


(36) 
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11,1 


(37) 


CRLB(Vb)  = 


From  the  last  equation  it  is  clearly  apparent  that  the  local  accuracy  in  the  bistatic  case  depends 
not  only  on  the  transmitted  waveform  but  also  on  the  bistatic  geometry  [Rih69],  [Tsa97]. 

It  is  important  to  observe  that  the  SNR  at  the  receiver  takes  into  account  the  energy  loss  due  to 
propagation: 

SNR  °c  — ]—  (38) 

krkt 


where  RT  =  J R ^  +  L2  +  2RsLsin  0R  is  the  range  from  transmitter  to  target. 

It  is  clear  that  for  L= 0  the  transmitter  and  the  receiver  are  co-located  and,  from  the  last 
equations,  it  is  clear  that  in  this  case  the  bistatic  FIM  coincides  with  the  monostatic  FIM. 

The  results  derived  in  this  Section  can  be  used  for  defining  a  tool  for  evaluate  the 
performance  of  a  given  monostatic  or  bistatic  channel  of  the  multistatic  system  and  for  design 
multistatic  weighting  coefficients  for  the  detection  process  that  will  be  described  in  the  next 
Section. 


6.2  Bistatic  Cramer- Rao  Lower  Bounds  for  a  burst  of  LFM  pulses 

Using  the  results  showed  in  previous  Section,  after  some  algebra,  it  is  possible  to  verify 
that,  in  the  monostatic  configuration,  the  FIM  of  a  burst  of  LFM  pulses  is  given  by  [Far09] 

fcVr2  k7r2T2 


iM{r,v)  =  -2SNR 


3 

k7t2T2 


jt2T2  7T2T2(\-N2) 


(39) 


Inverting  (39),  the  monostatic  CRLBs  for  the  delay  and  the  Doppler  are  given  by: 


CRLB(t)  : 


27t2T2k2SNR 


1  + 


f  rr  \ 


\TrJ 


n2-i 


(40) 


and 


CRLB(v) 


2n2T2SNR(N2  -\) 


(41) 
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These  results  are  in  agreement  with  those  obtained  in  [DogOl]. 

Using  this  result,  combined  with  (30),  it  is  possible  to  compare  the  monostatic  and  the  bistatic 
Root  Cramer-Rao  Lower  Bounds  (RCRLBs).  Figures  16-18  show  the  RCRLBs  of  range  and 
velocity,  both  in  the  monostatic  and  bistatic  case,  obtained  selecting  7=250pscc, 
Tr=  1msec,  B=  lMHz,/c=10GHz  and  N  =  8.  In  particular,  Figure  16  shows  the  RCRLBs  as  a 
function  of  the  receiver  to  target  range  Rr  when  Or  =0.  While  Figure  17  shows  the  results 
obtained  choosing  0r=  -0.4971.  Figure  18  shows  the  RCRLBs  as  a  function  of  the  receiver 
look  angle  Or  ,  both  in  the  case  of  Rr<L  and  Rr>L.  All  these  figures  have  been  obtained 
choosing  Vb=250  m/sec,  L=50km  and  holding  constant  the  SNR  to  OdB.  It  is  evident  that,  for 
all  the  parameter  values  we  tested,  the  bistatic  RCRLBs  are  always  higher  than  the  monostatic 
RCRLBs.  Anyway  when  the  distance  from  receiver  to  the  target  increases,  the  bistatic  system 
behaves  more  and  more  as  the  monostatic  one.  As  apparent  from  Figure  17,  the  effects  of 
geometry  are  prominent  where  the  target  approaches  the  baseline,  that  is  when  Rr<L  and  Or 
approaches  -n/2.  When  the  target  is  on  the  baseline,  the  RCRLBs  tend  to  infinity. 

In  this  case,  the  resulting  delay  is  Lie  and  the  radial  velocity  is  zero,  therefore  resolution  is 
totally  lost  and  the  RCRLBs  tend  to  infinity.  This  can  be  appreciated  by  realizing  that  the 
echo  arrives  at  the  receiver  at  the  same  instant  as  the  direct  signal,  independent  of  the  target 
location,  and  the  Doppler  shift  of  a  target  crossing  the  bistatic  baseline  must  be  zero,  because 
the  transmitter-to-target  range  changes  in  an  equal  and  opposite  way  to  the  target-to-receiver 
range,  independent  of  the  magnitude  and  direction  of  the  target  velocity.  However,  the  effects 
of  the  bistatic  geometry  are  less  prominent  when  the  distance  to  the  target  increases;  in  this 
case  the  bistatic  system  behaves  more  and  more  as  a  monostatic  system.  In  [TR4]  we  also 
calculated  in  closed  form  the  RCRLB  obtained  transmitting  a  Sinusoidal  Frequency 
Modulated  pulse.  This  signal  model  the  waveform  transmitted  by  a  FM  commercial  radio 
station  that  is  often  the  non  co-operative  transmitter  of  opportunity  exploited  by  a  Passive 
Coherent  Location  (PCL)  system.  The  obtained  results  are  pretty  similar  to  those  showed  in 
this  Section  and  are  not  reported  here  for  the  lack  of  space,  for  more  details  we  refer  the 
reader  to  [TR4] . 


-31  - 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


Figure  16  -  RCRLB  of  Range  and  Velocity  as  a  function  of  receiver  to  target  range  Rr, 
0r= 0;  L=50km,  SNR=0dB.  The  transmitted  signal  is  a  burst  of  LFM  pulses. 
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Figure  17  -  RCRLB  of  Range  and  Velocity  as  a  function  of  receiver  to  target  range  Rr, 
6r=  -0.49tt;  L=50km,  5VR=0dB.  The  transmitted  signal  is  a  burst  of  LFM  pulses. 


Figure  18  -  RCRLB  of  Range  and  Velocity  as  a  function  of  receiver  look  angle  0R; 
L=50km,  SNR=0dB.  The  transmitted  signal  is  a  burst  of  LFM  pulses. 
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7  Multistatic  Radar  Systems 

Multistatic  radars  utilize  multiple  transmitters  and  receivers.  Such  systems  differ  from 
typical  modem  active  radars  since  they  consist  of  several  different  monostatic  and  bistatic 
channels  of  observation.  Due  to  this  spatial  diversity,  these  systems  present  challenges  in 
managing  their  operation  as  well  as  in  usefully  combining  the  data  from  multiple  sources  of 
information  on  a  particular  area  of  surveillance.  The  information  gain,  obtained  through  this 
spatial  diversity,  combined  with  some  level  of  data  fusion,  can  give  rise  to  a  number  of 
advantages  over  both  the  individual  monostatic  and  bistatic  cases  for  typical  radar  functions, 
such  as  detection,  parameter  estimation,  tracking  and  identification.  As  showed  in  the 
previous  Sections,  the  performance  of  each  channel  of  the  multistatic  system  heavily  depends 
on  the  transmitted  waveform  and  on  the  geometry  of  the  scenario,  that  is,  the  position  of 
receivers  and  transmitters  with  respect  to  the  position  of  the  target. 

Exploiting  the  Monostatic  and  the  Bistatic  Cramer-Rao  Lower  Bounds,  it  is  possible  to 
calculate  the  channel  performance  of  each  TX-RX  pair  as  a  function  of  the  target  kinematic 
parameters.  In  particular,  in  this  Section  we  summarize  the  results  described  in  [TR4]  where 
we  presented  some  algorithms  to  select  the  best  channels  of  the  multistatic  system  and  to 
compute  the  rules  for  selecting  the  best  weighting  coefficients  for  fusing  the  signals  from 
multiple  receivers  in  order  to  improve  the  detection  performance  and  the  estimation  accuracy 
of  the  kinematic  parameters  of  the  target.  In  [TR4]  we  also  introduced  an  optimization 
methodology  for  selecting  only  some  channels  for  the  network,  independent  of  the  adopted 
fusion  rule.  All  the  techniques  described  in  this  Section  depend  on  the  Fisher  Information 
Matrix  (FIM)  which  is  itself  dependent  on  the  SNR,  the  AF  and  therefore  on  the  geometry  and 
the  transmitted  waveform. 


7.1  Optimal  channel  selection  in  a  multistatic  radar  system 

The  CRLB  study  carried  out  on  the  bistatic  geometry  can  be  applied  for  the  selection  of 
the  transmitter-receiver  (TX-RX)  pair  in  a  multistatic  radar  system.  We  have  seen  that  the 
performance  of  each  bistatic  channel  heavily  depends  upon  the  geometry  of  the  scenario  and 
the  position  of  the  target  with  respect  to  each  receiver  and  transmitter.  In  this  Section  we 
investigate  the  problem  of  optimally  selecting  the  TX-RX  pair,  based  on  the  information 
provided  by  the  CRLB  for  the  bistatic  geometry  of  each  bistatic  channel.  The  best  pair  is 
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defined  as  that  exhibiting  the  lowest  bistatic  CRLB  for  the  target  velocity  or  range  (or  a 
combination  of  the  two).  These  results  can  be  used  for  the  dynamical  selection  of  the  TX-RX 
signals  for  the  tracking  of  a  radar  target  moving  along  a  trajectory  in  a  multistatic  scenario. 

The  optimal  approach  for  a  multistatic  radar  system  would  be  central-level  track 
processing,  that  is,  to  send  all  observations  from  the  sensors  to  a  fusion  centre,  where  the 
observations  are  jointly  processed. 

This  approach  has  two  major  disadvantages: 

1)  the  amount  of  input  measurement  data  to  the  fusion  centre  may  be  very  high; 

2)  the  observations  from  the  sensors  with  the  worst  accuracy  and  resolution  can 
significantly  degrade  the  resolution  of  the  whole  systems. 

This  second  problem  is  maybe  the  most  important.  There  are  studies  that  show  that  the  best 
performance  of  the  fused  estimate  in  a  multisensor  system  occurs  when  the  sensors  have 
similar  accuracy  and  resolution  [Azi07] . 

If  accuracy  and  resolution  vary  widely,  there  is  a  risk  that  the  fused  track  performs  only 
marginally  better,  or  even  worse,  than  the  track  with  the  best  quality  estimate. 

As  an  example,  let  consider  the  case  of  a  target  which  is  close  to  the  baseline  joining  one  of 
the  transmitter  and  one  of  the  receiver  of  the  multistatic  system. 

In  this  case  the  range  and  Doppler  resolution  can  be  badly  degraded,  no  matter  what  the  radar 
waveform  is.  This  can  be  appreciated  by  realising  that  the  echo  arrives  at  the  receiver  at  the 
same  instant  as  the  direct  signal,  independent  of  the  target  location,  and  the  Doppler  shift  of  a 
target  crossing  the  bistatic  baseline  must  be  zero,  because  the  transmitter-to-target  range 
changes  in  an  equal  and  opposite  way  to  the  target-to-receiver  range,  independent  of  the 
magnitude  and  direction  of  the  target  velocity.  In  this  case,  resolution  is  totally  lost  and 
therefore  the  observation  from  this  transmitter-receiver  pair  could  hardly  degrade  the 
resolution  of  the  whole  multistatic  radar  system  [Azi07],  [GrelO],[Grell]. 

This  Section  proposes  an  optimization  algorithm  that,  in  a  generic  multistatic  scenario  and 
independent  of  the  adopted  fusion  technique,  specifies  what  channels  should  be  discarded  and 
what  channels  should  be  considered  during  the  fusion  process. 

Using  this  algorithm,  only  a  subset  of  data  are  communicated  to  the  fusion  centre,  more 
specifically  only  from  those  sensors  exhibiting  the  best  performance  in  terms  of  estimation 
accuracy  of  the  target  parameters. 

In  the  most  general  case,  the  multistatic  scenario  is  the  one  pictorially  depicted  in  Figure  19, 
where  there  are  M  transmitter  and  N  receivers,  co-located  or  not,  surveying  a  common 
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coverage  area.  It  is  supposed  that  a  set  of  orthogonal  waveform  is  transmitted,  where 
orthogonality  is  assumed  to  be  maintained  for  any  Doppler-delay  shift.  With  proper  design, 
transmit-receive  paths  can  be  separated,  so  that  each  sensor  can  receive  echoes  due  the  signals 
generated  by  all  the  transmitters  and  can  select  the  transmitted  signal  of  interest.  In  this  way, 
the  multistatic  system  can  be  considered  as  consisting  of  NM  different  monostatic  and  bistatic 
channels  of  observation.  The  proposed  optimization  algorithm  exploits  the  results  obtained  in 
the  previous  Sections  to  approach  the  problem  of  optimally  selecting  the  channels  to  be  used 
by  the  fusion  process.  This  algorithm  is  divided  in  two  steps.  In  the  first  step,  each  of  the  N 
receivers  of  the  network  selects  one  of  the  M  transmitters  in  order  to  obtain  the  best 
performances  in  terms  of  estimation  accuracy  of  the  target  range  and/or  velocity.  While,  in  the 
second  step,  after  ranking  the  so  obtained  N  channels  from  the  worst  to  the  best,  where  the 
best  is  the  one  that  exhibits  the  lowest  bistatic  CRLB,  only  the  first  N'<N  are  selected  for  the 
fusion  process.  In  the  following,  with  the  help  of  an  illustrative  example,  we  describe  in  detail 
how  the  optimization  algorithm  works.  In  the  first  step  of  the  optimization  process,  each 
receiver  selects  the  best  transmitter  on  its  own.  In  a  multistatic  network  consisting  of  M 
transmitters  and  N  receivers  for  each  receiver  there  are  M  different  channels.  In  the  case  that 
the  considered  receiver  is  colocated  with  one  of  the  transmitter,  there  are  M- 1  bistatic 
channels  and  1  monostatic  channel. 


Tx  3  RX  N 

•  ♦ 


Rx  3 

♦ 


O  Target 
0  Transmitter 
♦  Receiver 


Figure  19  -  Generic  Multistatic  Scenario. 
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As  an  example,  let  consider  a  surveillance  map  of  dimension  Lx=20km  and  Ly=20km  where 
there  are  seven  transmitters  placed  at  coordinates  T1 1  '=(1 3.9  km, 15.9  km), 
T(2)=(6.3  km,  3.7  km),  T(3)=(19  km, 9.8  km),  T(4)=(0.7  km, 8.9  km),  T(5,=(8.8  km,  12.9  km), 
T<6)=(7.6  km,  14.2  km),  T(7,=(15.3  km,  15.1  km)  and  the  receiver  is  co-located  with  the  first 
transmitter.  In  this  example  it  is  supposed  that  the  transmitters  send  orthogonal  signals  with 
the  same  power  and  characterized  by  the  same  AF.  Therefore,  the  different  performances 
among  the  channels  depend  only  on  the  propagation  path  loss  and  the  configuration  of  the 
bistatic  geometry.  For  each  point  of  the  analysed  area,  evaluating  the  performance  of  each 
bistatic  channel,  it  is  possible  to  select  the  transmitter  with  the  best  performance,  that  is,  the 
one  with  the  lowest  CRLB  on  the  target  range  and  velocity  estimation  accuracy. 

Figures  20  (a)  and  (b)  show,  in  a  colour  coded  map,  the  transmitter  to  be  selected  in  order  to 
provide  the  minimum  CRLB  for  each  point  of  the  analysed  area  for  range  and  velocity 
estimation,  respectively.  The  color-map  of  these  figures  is  clustered  into  7  colours,  each  of 
which  is  associated  with  one  of  the  7  transmitters.  It  is  apparent,  that  the  results  are  very 
similar,  but  it  is  also  possible  to  build  a  cost  function  using  a  weighted  combination  of  the  two 
CRLB.  Figures  21  (a)  and  (b)  show  the  minimum  Root  of  the  CRLB  (RCRLB)  of  the  target 
range  and  the  target  velocity,  respectively,  that  is  the  value  of  the  RCRLB  which  is  provided 
by  the  transmitter-receiver  pair  which  has  the  minimum  RCRLB. 

From  this  example,  it  is  apparent  that  each  receiver  knows,  for  each  point  of  the 
surveillance  map,  which  is  the  transmitter  to  be  selected  in  order  to  meet  the  best 
performance.  Therefore,  as  pictorially  shown  in  Figure  22,  based  upon  the  actual  estimate  of 
the  target  position,  each  receiver  can  dynamically  select  the  signal  of  the  best  transmitter  and 
discard  the  signals  transmitted  by  the  other  sensors.  Doing  so,  only  N  of  the  MN  possible 
channels  are  selected  for  the  fusion  process. 

Even  so,  from  Figures  21  (a)  and  (b)  it  is  apparent  that,  for  same  points  of  the  surveillance 
map,  the  best  performance  that  a  receiver  can  achieve  could  be  poor.  This  is  only  due  to  the 
geographical  distribution  of  the  sensors  in  the  multistatic  system  but,  for  the  reasons 
previously  described,  this  could  be  a  problem  for  the  fusion  process. 

This  problem  can  be  solved  by  selecting  a  subset  of  channels  among  the  N  previously 
obtained.  In  particular,  in  the  second  step  of  the  optimization  algorithm,  the  N  channels  are 
ranked  from  the  worst  to  the  best  and  only  the  first  N’  are  selected  for  the  fusion  process.  The 
number  N’  could  be  also  dynamically  changed  by  the  fusion  process  in  order  to  meet  pre¬ 
specified  performance  goals. 
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Figure  20  -  (a)  Optimum  Transmitter  for  target  range  estimation;  (b)  Optimum 
Transmitter  for  target  velocity  estimation.  The  receiver  is  co-located  with  T(1). 


Figure  21  -  (a)  Minimum  RCRLB  of  the  target  range  [dBm];  (b)  Minimum  RCRLB  of 

the  target  velocity  estimation  [dBm/sec]. 


Tx  1 


Figure  22  -  Optimum  transmitter  selection. 
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7.2  Optimum  weighting  rule  for  multistatic  detection 

This  Section  deals  with  the  description  of  a  coherent  multistatic  system  where  the  rule  for 
selecting  the  weighting  coefficients  is  computed  exploiting  the  results  obtained  in  previous 
Sections.  For  ease  of  notation,  the  described  system  is  composed  of  one  transmitter  and  N 
receivers.  In  the  case  of  multiple  transmitters,  it  is  supposed  that  it  is  possible  to  separate  each 
transmitter-receiver  path  so  that  each  receiver  can  select  the  transmitted  signal  of  interest.  In 
this  way,  the  multistatic  system  is  composed  by  MN  channels  of  observation,  where  M  is  the 
number  of  transmitters.  The  only  difference  with  the  case  of  a  single  transmitter  is  only  on  the 
number  of  channels,  but  the  algorithm  remains  the  same.  The  multistatic  receiver  that  will  be 
implemented  refers  to  a  centralized  nomenclature  where  the  data  at  each  receiver  are  first 
collected  by  a  central  processor  and  then  jointly  combined  for  the  detection. 

In  particular,  the  central  processing  of  the  analyzed  system  will  exploit  part  of  the  results 
described  in  Section  7.1.  In  particular,  the  optimization  algorithm  described  in  Section  7.1  can 
be  viewed  as  a  fusion  algorithm  that  uses  hard  weighting  coefficients,  that  is,  weights  that  can 
assume  only  two  values:  zero  in  the  case  of  a  bad  channel  and  1  in  the  other  case.  On  the  other 
hand,  the  weighting  coefficients  of  the  central  processor  described  in  this  Section  are  soft,  that 
is  they  can  assume  values  that  belongs  to  the  range  [0,1]  depending  on  the  information  gained 
by  the  Fisher  Information  Matrix. 

Figures  23  and  24  show  the  block  diagram  of  the  multistatic  receiver  that  will  be 
described  in  this  Section.  From  Figure  23,  it  is  apparent  how  the  detection  processing  is 
divided  in  two  parts:  the  local  processing  at  each  receiver  and  the  central  processing  at  the 
central  processor.  Figure  24  shows  the  flow  diagram  of  the  signal  processing  at  the  central 
processor. 

As  showed  in  Figure  23,  the  flow  diagram  of  each  receiver  of  the  multistatic  network  is 
the  most  adopted  in  modem  coherent  radar  system  where,  after  the  classical  processing  to 
down-shift  the  received  signal  to  the  base  band,  the  receiver  performs  a  range  -  Doppler 
matched  filtering.  In  the  case  of  multiple  transmitters,  the  range  -  Doppler  matched  filtering 
consists  in  a  bank  of  filters,  each  of  which  is  matched  to  the  waveform  emitted  by  each 
transmitter.  As  previously  discussed,  the  performance  of  each  receiver  depends  on  the 
transmitted  waveform  but  unfortunately  is  also  heavily  sensitive  to  the  geometry,  that  is,  the 
position  of  the  receiver  and  the  transmitter  with  respect  to  the  position  of  the  target. 

For  a  fixed  geometry  of  the  simulated  scenario,  it  is  possible  to  derive  easily  the  CRLBs 
of  each  channel  and  therefore  its  performance. 
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Figure  23  -  Flow  Diagram  of  the  multistatic  detector  [TR4]. 
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Figure  24  -  Flow  Diagram  of  the  Central  Processor. 


Moreover,  it  is  clear  that,  each  receiver  of  the  multistatic  system  has  its  own  local 
geometry,  that  is,  each  receiver  knows  at  what  range  from  its  location  there  is  a  target  and  the 
target  speed  straight  towards  or  away  from  the  receiver. 

The  role  of  the  central  processor  is  therefore  to  exploit  the  information  gained  by  the 
knowledge  of  the  geometry  of  the  network  (the  position  of  each  transmitter  and  each  receiver) 
to  convert  the  local  coordinates  of  the  receivers  into  global  coordinates. 
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The  information  about  range  and  radial  velocity  from  each  receiver  are  therefore  converted 
into  information  about  the  target  state,  that  is,  the  real  position  of  the  target  and  the  velocity 
components  along  the  axes  of  the  reference  coordinates  system.  Moreover,  another  important 
role  of  the  central  processor  is  to  evaluate  for  each  point  of  the  global  surveillance  map,  that  is 
for  each  cell  (x,y,Vx,Vy),  the  performance  that  can  be  gained  for  each  path  of  the  multistatic 
system.  Therefore  the  role  of  the  central  processor  is  to  generate  a  global  coordinates 
detection  map  exploiting  the  information  about  the  geometry  of  the  multistatic  network. 

Let’s  describe  now  each  block  of  the  flow  diagram  of  Figure  24.  The  block  buffer  indicates 
that,  for  each  scan,  the  central  processor  stores  the  local  range-Doppler  maps  of  each  receiver. 
In  particular,  each  local  map  can  be  viewed  as  a  two  dimensional  matrix  whose  elements  are 
the  samples  at  the  output  of  the  matched  filter  of  each  receiver. 

The  central  processor  scans  each  point  of  the  global  surveillance  area  described  by  the  vector 
(x,y,Vx,Vy),  that  is  the  Cell  Under  Test  (CUT).  Once  fixed  the  CUT  and  knowing  the 
geometry  of  the  multistatic  system,  the  role  of  the  block  range-Doppler  selection  is  to  selects 
the  corresponding  range-Doppler  cell  of  each  receiver. 

Thus,  the  problem  of  detecting  a  target  signal  in  a  fixed  resolution  cell  of  the  surveillance 
map  can  be  posed  in  terms  of  the  following  binary  hypotheses  test: 


where  the  sample  r;-  is  the  output  of  the  filter  matched  to  the  normalized  transmitted  signal  for 
the  i'-th  receiver  in  the  range-velocity  cell  corresponding  to  the  point  (x,y,Vx,Vy).  Note  that  if 
the  transmitted  signal  is  a  burst  of  pulses,  the  sample  r;-  is  the  output  of  the  filter  matched  to 
the  entire  burst.  In  the  signal  model  in  eq.  (42),  n,  is  the  additive  noise  at  the  /-th  sensor;  it  is 
modelled  as  a  zero-mean  complex  Gaussian  variable,  in  short  n,  ~  CN( 0,  a2).  The  samples  n, 
are  assumed  to  be  independent  from  channel  to  channel  and  identically  distributed.  The 
parameter  a,  accounts  for  the  channel  propagation  effect  and  the  target  radar  cross  Section. 
We  considered  the  case  in  which  the  amplitude  a,  is  a  random  variable.  It  depends  on  the 
target  bistatic  RCS  of  the  /-th  bistatic  angle.  This,  together  with  the  mentioned  features  of  the 
multistatic  geometry,  justifies  the  assumption  that  the  ay  are  mutually  independent  random 
quantities.  With  regard  to  the  marginal  PDFs,  a ;  is  a  zero-mean  complex  Gaussian  variable 
whit  variance  a 2  varying  from  path  to  path  a,-  ~  CN(0,  a,2).  In  particular  the  Signal  to  Noise 
power  Ratio  at  the  input  of  the  i- th  channel  is  defined  as  SNRi=  a2/  o2  =  Const/(RT{,)RR{’>j2  and, 
as  previously  mentioned,  is  inversely  proportional  to  the  energy  path  loss  factor  ( Rtu>Rr(i) )2 
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due  to  propagation  where,  Rt']  and  RrU)  are  the  range  from  transmitter  to  target  and  the  range 
from  receiver  to  target  for  the  /-th  channel,  respectively.  The  parameter  c;-.  that  appears  in  eq. 
(42)  is  a  complex  number  whit  absolute  value  lower  than  or  equal  to  one  which  accounts  for 
the  other  effects  of  propagation  and  scattering  along  the  /-th  path.  The  meaning  of  c,  will  be 
discussed  later  on,  for  the  moment  let  consider  the  case  in  which  c,=l.  Under  this  assumption 
and  considering  that  the  amplitude  a-,  and  the  noise  n,  are  independent,  the  observation  r;  is  a 
zero-mean  complex  Gaussian  variable  under  both  hypotheses.  Therefore,  it  is  possible  to 
write: 

rt  ~CN(Q,<t?  +a2)  Hl 

,  '  (43) 

r  -  C/V(0,  (T2)  H0 


Hence,  by  exploiting  the  independence  of  individual  likelihoods  in  each  channel,  it  is  possible 
to  verify  that  the  Neyman-Pearson  decision  test  develops  into  the  following  form 
[Dad86],[Con83]: 

Eft H  $2  (44) 

i= 1  «i 


where 

SNR, 

Pl  ~  1  +  SNR, 


(45) 


Note  that  the  weight  /?,  are  non  negative  and  are  an  increasing  function  of  SNRt.  therefore  the 
central  processor  emphasizes  those  channels  along  which  the  SNRs  are  the  highest.  Note  that 
the  receiver  to  implement  the  given  test,  depending  on  cell  under  test  of  the  surveillance  map, 
needs  to  continually  update  the  weights,  which  are  themselves  dependent,  through  SNRj,  on 
the  distances  from  the  transmitter  to  the  target  and  from  the  target  to  the  receiver. 

Consider  now  the  case  in  which,  under  the  H\  hypothesis,  c;  in  eq.  (42)  is  a  complex 
random  variable.  This  simple  assumption  has  been  done  in  order  to  model  the  effects  of  the 
bistatic  geometry  along  the  /-th  path.  In  particular  we  generated  c,  as 

ci=X(Ti,Vi)  (46) 


where  r,=  zHr  U  and  v;-  =  vm  -  va  are  the  delay  and  Doppler  shift  for  the  z-th  path  obtained 
using  equations  (23)  and  (24)  and  by  random  generating  the  estimates  of  the  range  and 
velocity  as 

R$)~N(R%\CRLBf)  (47) 
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Vg(0  ~  N{Vgn  ,  CRLB^ )  (48) 

where  and  VRn  are  the  actual  range  and  bistatic  velocity  of  the  target,  while  CRLB^’  and 
CRLBj!1  are  the  Cramer-Rao  lower  bounds  of  range  and  velocity  in  the  i-th  path. 

By  random  generating  and  VRl) ,  the  value  of  c,  takes  into  account  the  miss-matching  at 

the  i-th  receiver,  moreover,  due  the  geometry  dependent  non  linear  transformation  of 
equations  (23)  and  (24),  that  gives  the  delay  and  the  Doppler  as  a  function  of  the  range  and 
the  bistatic  velocity,  the  value  of  c,  takes  also  into  account  the  effects  of  the  distortion  of  the 
Ambiguity  Function  due  to  the  bistatic  geometry. 

As  an  example,  let  consider  the  case  in  which  CRLB','1  and  CRLBj!'  are  low,  in  this  case  the 

values  assumed  by  R{r  and  VR  ‘  are  near  the  actual  range  and  velocity,  therefore  the  values 
assumed  by  r,  and  v,  are  almost  zero,  and  hence  the  value  assumed  by  c,  is  near  the  maximum 
that  is  1.  On  the  other  hand,  for  a  bad  bistatic  geometry,  the  values  assumed  by  CRLB^’  and 

CRLBj:'  are  very  high,  therefore  there  is  a  high  probability  that  also  r,  and  v,  are  high  and 

hence  the  value  assumed  by  them  is  almost  zero  and  the  observed  signal  is  only  noise. 

In  other  words,  by  generating  in  this  way  the  value  of  c„  the  signal  model  in  eq  (42)  is  now 
dependent  on  the  sharpness  of  the  CAF  around  its  maximum  and  to  the  distortion  of  its 
behaviour  due  to  the  bistatic  geometry.  It  is  clearly  apparent  that  it  is  very  difficult  to  derive 
the  PDF  of  c,  and  hence  it  is  difficult  to  derive  the  Neyman-Pearson  decision  test.  By  the  way, 
our  approach  is  to  use  the  same  receiver  of  eq.  (44)  but  choosing  the  weights  pt  in  a  different 
manner.  In  fact,  the  weights  in  eq.  (45)  depends  only  on  the  energy  path  loss  and  they  do  not 
take  into  account  the  other  effects  due  to  the  bistatic  geometry.  As  an  example,  let  consider 
the  case  of  a  target  is  in  the  baseline  between  the  transmitter  and  the  i-th  receiver.  When  the 
distance  between  the  transmitter  and  i-th  receiver  is  low,  the  expected  signal  to  noise  ratio  is 
high  and  therefore  the  central  processor  using  the  weights  in  (45)  tends  to  emphasize  this  path 
during  the  detection  process.  As  showed,  in  this  case  the  resolution  is  totally  lost  and 
therefore  the  observations  from  this  path  can  significantly  degrade  the  performance  of  the 
whole  system.  Therefore,  the  rule  for  selecting  the  weighting  coefficients  should  be  different. 
In  particular,  the  weights  should  be  highlight  those  channels  that  exhibit  the  best  performance 
in  terms  of  estimation  accuracy  of  the  target  parameters  instead  of  emphasizes  those  along 
which  the  SNRs  are  the  highest.  In  this  work  we  propose  the  following  weighting  rule 


-43  - 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


(49) 


P  ^  - 

r+Tmce{yB'(R"'.V«’)} 

where  Trace {.}  is  the  trace  operator  and  y  is  a  constant  that  we  fixed  equal  to  the  inverse  of 
the  sum  of  the  CRLBs  of  range  and  velocity  in  the  monostatic  case  when  SNR =0dB.  In  this 
way  the  weighting  coefficients  depend  on  the  SNR  and  the  normalization  constraint  can  be 
viewed  as  a  terms  of  comparison  of  the  performance  of  the  bistatic  channel.  That  is,  a  channel 
can  be  considered  a  good  channel  if  trace  {Jb"1}  is  lower  than  y'  ,  while  a  channel  is  bad  if 
trace  {Jg  1 }  is  greater  than  y"1.  It  is  important  to  recall  that  R^'  and  VB(,)  are  the  actual  range 
and  velocity  of  the  target.  Therefore,  depending  on  the  cell  (x,y,Vx,Vy)  under  test,  the  central 
processor  can  easily  calculate  these  two  values  for  each  channel  knowing  the  position  of  each 
sensor  of  the  multistatic  system. 

As  apparent  in  eq.  (49),  the  weights  p,  depend  on  the  inverse  of  the  bistatic  Fisher  Information 
Matrix  which  is  dependent  on  the  energy  path  loss  through  SNR,,  on  the  geometry  and  also  on 
the  transmitted  waveform.  This  is  also  important  in  the  case  of  a  multistatic  system  where 
multiple  transmitters  are  used.  Using  the  weights  in  (49)  the  central  processor  emphasizes 
those  channels  that  exhibit  the  best  performance  in  terms  of  estimation  accuracy  and  is  also 
able  to  discard  those  channels  where  the  resolution  is  totally  lost.  In  this  case  the  trace  of  the 
inverse  of  the  bistatic  FIM  tends  to  infinity  and  the  therefore  the  weight  p ,•  tends  to  zero. 

Note  that  the  value  of  y  can  be  fixed  to  the  inverse  of  the  sum  of  the  CRLBs  of  range  and 
velocity  in  the  monostatic  case  when  SNR=0dB.  In  this  way  the  weighting  coefficients  depend 
on  the  SNR  and  the  normalization  constraint  can  be  viewed  as  a  terms  of  comparison  of  the 
performance  of  the  bistatic  channel.  That  is,  a  channel  can  be  considered  a  good  channel  if 
trace  {Jg1 }  is  lower  than  y'1,  while  a  channel  is  bad  if  trace  {Jg1 }  is  greater  than  y'1. 

Note  that  the  receiver  to  implement  the  given  test  for  both  the  weighting  rules,  depending  on 
cell  under  test  of  the  surveillance  map,  needs  to  continually  update  the  weights.  Moreover, 
note  that  the  matrix  in  eq.  (49)  is  2x2  and  therefore  the  inversion  is  straightforward. 
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7.3  Performance  of  the  optimum  multistatic  detector 

Figure  25  shows  the  multistatic  scenario  analyzed  in  our  simulations.  The  multistatic 
system  is  composed  by  one  transmitter  and  three  receivers.  The  distance  between  the 
transmitter  and  each  receiver  is  the  same  and  equal  to  30km.  The  transmitted  signal  is  a  burst 
of  Linear  Frequency  Modulated  pulses  where  T=250psec,  Tr=  1msec,  B=  1MHz, /c=10GHz 
and  A  =  8  (Number  of  coherent  pulses). 

As  showed  in  Figure  25,  in  our  simulation  we  considered  three  cases.  In  the  first  one  we 
considered  a  target  in  the  middle  of  the  first  baseline,  that  is  between  the  transmitter  and  the 
first  receiver.  In  this  case  the  target  is  moving  with  a  speed  of  250  m/sec  and  direction  of  135° 
with  respect  to  the  horizontal  axis.  In  the  second  case  the  target  is  in  the  second  baseline,  10 
km  far  from  the  transmitter,  moving  with  a  speed  of  250  m/sec  and  direction  of  150°.  While 
in  the  third  case,  the  target  is  15  km  far  from  the  third  baseline  and  20  km  far  from  the  second 
baseline,  and  it  is  moving  with  a  speed  of  250  m/sec  and  direction  of  180°  whit  respect  to  the 
horizontal  axis.  Table  1  shows  the  weighting  coefficients  obtained  with  both  the  methods,  that 
is  the  “ SNR' ’  rule  in  eq.  (45)  and  the  ‘Trace”  rule  in  eq.  (49).  Note  that  the  coefficients  have 
been  normalized  in  order  to  satisfy  the  relationship  Z/?,=  1 . 

As  apparent  from  Figure  25,  in  Case  1  the  target  is  in  the  first  baseline,  therefore  the 
resolution  of  the  first  receiver  is  totally  lost  also  if  the  SNR  for  this  path  is  the  highest. 

As  apparent  from  Table  1,  the  coefficient  corresponding  to  the  first  receiver  is  null  using  the 
Trace  rule  while  it  is  the  highest  using  the  SNR  rule.  Similar  considerations  can  be  drawn  for 
Case  2.  In  the  third  case  the  target  is  in  an  optimal  position  for  all  the  receiver,  the  main 
difference  from  one  receiver  to  the  other  is  related  only  to  the  energy  path  loss,  therefore  the 
weighting  coefficients  obtained  with  both  the  rule  are  almost  the  same.  Figure  26  shows  the 
Receiver  Operating  Characteristic  (ROC)  of  the  multistatic  detector  in  (44)  for  all  the  three 
analyzed  cases.  We  also  plotted  the  ROC  obtained  by  forcing  all  the  coefficients  to  1/A. 

As  apparent  from  the  results  the  ROC  obtained  in  Case  3  (dash-dotted  lines)  are  almost  the 
same  for  the  three  weighting  rules,  while  in  Case  1  (solid  lines)  and  Case  2  (dashed-lines)  the 
performances  obtained  whit  the  Trace  rule  are  the  best.  Moreover,  in  Case  1,  the  performance 
obtained  with  the  SNR  rule  is  lower  than  that  obtained  with  all  the  weights  equal  to  1/A.  In 
this  case,  the  SNR  for  the  first  receiver  is  the  highest  but,  due  to  bad  geometry,  the  resolution 
is  totally  lost.  Using  the  SNR  rule,  the  central  processor  emphasizes  the  observations  from  the 
first  sensors  and  this  strongly  degrades  the  performance  of  the  multistatic  detector. 
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Probability  of  Detection 


SNR 

TRACE 

Case  1 

p=[0.3495;  0.3328;  0.3177] 

p= [0.0000;  0.4984;  0.5016] 

Case  2 

p=[0.3317;  0.3366;  0.3317] 

p=[0.5000;  0.0000;  0.5000] 

Case  3 

p=[0.2735;  0.3553;  0.3712] 

p=[0.2939;  0.3449;  0.3612] 

Table  1-  Weigting  Coefficients. 
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Figure  25  -  Simulated  scenario. 
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Figure  26  -  Receiver  Operating  Characteristic  of  the  Multistatic  Detector. 
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8  Passive  Coherent  Location  (PCL)  Systems 

Recently,  great  interest  has  been  devoted  to  Passive  Coherent  Location  (PCL)  systems, 
that  is,  bistatic  radars  that  exploit  the  signal  emitted  by  illuminators  of  opportunity,  such  as 
broadcast  or  communications  signals.  PCL  systems  have  some  significant  attractions,  in 
addition  to  those  common  to  all  bistatic  radars:  the  PCL  receivers  do  not  need  any  transmitter 
hardware  of  their  own,  are  completely  passive,  and  hence  undetectable.  Among  all  the 
potential  signals  of  opportunity  we  can  mention  FM  radio  and  Digital  Audio  Broadcast 
(DAB)  [DiL08],  [Mal07],  [How05],  [GrelO],  analogue  TV  and  Digital  Video  Broadcast 
(DVB)  [And05]  [Gri86],  [How99],  Wi-Fi  [Guo08]  or  cell  phone  network  signals  [Stil2], 
[Sun08],  [DeM08],  [Men06],  [Pet09], 

Cellular  phone  networks  are  now  ubiquitous  in  most  countries  and  the  implementation 
and  application  of  PCL  systems  that  exploit  these  signals  is  an  active  research  area.  Over  the 
last  few  years  there  has  been  an  evolution  of  mobile  networks  towards  the  third  generation 
radio  wireless  communications  (3G)  such  as  UMTS  (Universal  Mobile  Telecommunications 
System).  The  UMTS  is  the  European  standard  for  3G  wireless  communication  and  it  has  been 
developed  by  3GPP  (3rd  Generation  Partnership  Project)  [TS06].  The  downlink  signal 
transmitted  from  the  UMTS  base-station  presents  interesting  properties  that  suggest  the 
usability  of  the  UMTS  signal  as  opportunity  signal  for  passive  radar  systems.  These  properties 
are:  the  large  bandwidth,  that  implies  a  high  range  resolution;  the  use  of  orthogonal  codes, 
that  guarantees  low  sidelobes  of  the  Ambiguity  Function;  and  the  growing  coverage  of  UMTS 
signals  on  the  international  territory,  which  renders  a  multistatic  radar  configuration  feasible. 
Moreover,  the  UMTS  system  uses  a  band  that  spans  from  2110-2170  MHz,  that  makes  easy  to 
work  with  directive  antennas. 

In  particular,  the  UMTS  signal  has  the  following  characteristics  [TS06],  [Wal98], 
[Wal03].  There  are  two  forms,  Frequency  Division  Duplex  (FDD)  and  Time  Division  Duplex 
(TDD),  each  of  which  has  a  bandwidth  of  5  MHz.  FDD  requires  two  frequency  bands,  one  for 
the  up-link  and  one  for  the  down  link,  while  the  TDD  requires  a  single  band.  The  transmission 
is  Wideband  CDMA  (WCDMA)  using  Walsh-Hadamard  coding  with  a  transmission  rate  of 
3.84  Mchips/s.  The  modulation  used  is  QPSK  (Quadrature  Phase  Shift  Keying)  where  the 
pulses  are  shaped  with  a  Root  Raised  Cosine  (RRC)  filter  with  roll-off  factor  of  0.22.  The 
study  of  PCL  systems  that  exploits  the  downlink  signal  of  a  UMTS  base  station  is  of  great 
interest,  anyway  the  performance  of  such  systems  has  been  evaluated  only  with  some 
experimental  systems  [Pet09],  [Gri03],  [Gri05]. 
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In  this  Section  we  evaluate  the  performance  of  these  systems  by  calculating  the 
Ambiguity  Function  (AF)  of  a  QPSK  signal  and  the  modified  Cramer- Rao  lower  bounds 
(MCRLBs)  on  the  estimation  accuracy  of  target  range  and  velocity,  for  both  monostatic  and 
bistatic  configurations.  These  results  provide  a  useful  tool  for  the  evaluation  of  the 
performance  of  a  PCL  system  that  exploits  the  downlink  signal  of  a  UMTS  base  station  and 
they  can  also  be  used  to  design  an  algorithm  for  the  dynamical  selection  of  the  transmitter  in  a 
multistatic  network  where  the  receiver  can  exploit  the  signals  emitted  from  several  base 
stations. 


8.1  The  UMTS  downlink  transmitted  signal 

The  downlink  signal  transmitted  by  a  UMTS  base  station  is  a  QPSK  signal  where  the 
pulses  are  shaped  with  a  RRC  filter  [Wal98],  [Wal03],  [TS06]. 

The  unit  energy  complex  envelope  of  the  QPSK  signal  is  given  by: 

i  N- 1 

u(t)  =  -j=Yjcng(t-nT),  (50) 

\N  n= 0 


where  c„  are  the  transmitted  symbols,  N  is  the  number  of  transmitted  symbols,  and  T  is  the 
inverse  of  the  symbol  rate  (or  chip  rate).  In  the  QPSK  modulation,  symbols  c„  belong  to  the 
set  {  e  7‘2''+l,ir/4,  /z=  ]  ,2,3,4 }  and  they  are  assumed  to  be  independent  and  identically  distributed 

(HD): 


n  =  k, 
n  ^  k. 


(51) 


The  pulse  git)  is  a  delayed  RRC  pulse,  in  particular: 


git)  =  h 
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(52) 


where  the  delay  D  is  chosen  such  that 
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is  a  centered  RRC.  The  parameter  a  is  the  roll-off  factor,  with  0<a<l.  The  Fourier  Transform 
(FT)  of  the  pulse  h(t)  is  equal  to: 

H{f)  =  4ciT)  (55) 
where  C(f)  is  the  FT  of  a  Raised  Cosine  (RC)  pulse  c(t).  C(f )  and  c{t )  are  defined  as: 

.  1  -a 
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c(t)  =  IFT  {C(/)}  =  sine 
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(57) 


where  IFT{-}  is  the  Inverse  Fourier  Transform  (IFT)  operator.  In  Figs.  27  and  28,  hit)  and 
H(f)  are  plotted  for  different  values  of  the  roll-off  factor  a.  It  is  easy  to  verify  that: 


\h2(t)dt=\H2(f)df=\C(f)df  =  1 


(58) 


In  the  following  sections  we  evaluate  the  AF  and  the  MCRLBs  of  the  QPSK  signal.  These 
quantities  are  then  used  to  measure  the  performance  of  a  PCL  system  that  exploits  a  UMTS 
base  station  as  a  transmitter  of  opportunity.  In  particular,  for  this  kind  of  system  the  value  of 
the  roll-off  factor  is  equal  to  0.22  while  the  inverse  of  the  symbol  rate  T  is  equal  to  0.26  ps 
[Wal98],  [Wal03],  [TS06], 

In  order  to  achieve  a  fine  Doppler  resolution,  PCL  radar  systems  have  long  observation 
time,  typically  of  the  order  of  0.1  s.  Hereafter,  without  loss  of  generality,  we  will  fix  the 
observation  time  NT  to  0.1  s,  in  this  case  the  corresponding  number  of  integrated  pulses  is 
A=384615. 
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Figure  27  -  Impulse  response  of  a  RRC  for  different  values  of  the  roll-off  factor  a. 


Figure  28  -  Frequency  response  of  a  RRC  for  different  values  of  the  roll-off  factor  a. 
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8.2  Monostatic  and  Bistatic  Ambiguity  Function  of  the  UMTS  signal 

Referring  to  the  definition  of  CAF  in  eq.  (1),  and  reported  in  (59)  for  convenience,  we 
have  that,  for  a  fixed  sequence  of  symbols  {c„,  //=(),  1 , . . .  ,N- 1  ( .  the  CAF  of  the  transmitted 
pulse  in  (1)  is  defined  by  [Tsa97]: 

+oo 

ZAT>V)=  J  u(t)u* (t-r)e~j2m'dt ,  (59) 


where  t=th-t0  is  the  difference  between  the  hypothesized  delay  and  the  actual  delay,  while 
v=vH-va  is  the  difference  between  the  hypothesized  Doppler  shift  and  the  actual  Doppler  shift. 
The  CAF  in  (59)  describes  the  properties  of  the  signal  for  a  fixed  sequence  of  symbols.  To 
evaluate  the  target  resolution  of  the  transmitted  waveform  independently  of  the  particular 
sequence  of  transmitted  symbols,  we  can  evaluate  the  average  CAF  (ACAF)  defined  as: 

x(r,v)  =  E{xc(z,v)},  (60) 


where  the  expectation  £{•)  is  with  respect  to  the  vector  of  symbols  c  =[c i,  ...,  cw]. 
By  these  considerations,  the  CAF  of  the  transmitted  signal  is  given  by 
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(61) 


where  we  used  the  IID  assumption  for  the  symbols  in  eq.  (51)  and  /g(r,v)  is  the  CAF  of  g(t) 
given  by: 

Xg(r,v)  =  ej*Dvx„(T,v),  (62) 

with 

+oo 

Xh  (t,v)  =  J  h(t)h* (t -  T)e~j27rvtdt  =  FT{h{t)h\t-  z)}  =  H(v)®  H(v)ej2™  (63) 
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where  ®  is  the  convolution  along  the  v  axis  while  FT{  • }  the  Fourier  Transform  operator.  The 
average  AF  (AAF)  is  defined  as  the  absolute  value  of  the  average  CAF.  Hence,  the  AAF  of 
the  transmitted  pulse  is  given  by: 


z(av) 


1  sin(;r./V7v) 
N  sin  (itTv) 


Z,(Av) 


(64) 


and  exhibits  the  maximum  in  (t,  v )  =  (0,  0).  The  AAF  is  plotted  in  Fig.  29  for  NT=  0.1  s, 
T=  0.26  ps,  and  a=0.22. 


|X(T,V)| 


Figure  29  -  AAF  of  a  QPSK  signal  with  NT=0.1  s,  T=0.26  ps,  and  a=0.22. 
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To  evaluate  the  delay  and  Doppler  resolution  of  the  transmitted  pulse,  the  zero-Doppler 
and  zero-delay  cuts  of  the  AAF  are  often  analyzed.  The  zero-Doppler  cut  is  the  AAF  at  v=0, 
that  is  the  autocorrelation  function  of  the  transmitted  pulse.  It  represents  the  matched  filter 
output  when  there  is  no  Doppler  mismatch.  Using  (64)  and  the  definition  of  the  RRC  pulse,  it 
is  easy  to  verify  that: 


\z{t,0)\  =  | %h  (t,0)|  =  | h(T) ®  h(T) I  =  \lFT {C(v)}| 


sine 


cos 

'  a 7tt\ 
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On  the  other  hand,  the  zero-delay  cut  is  the  AAF  for  r  =0  and  represents  the  output  of  the 
matched  filter  at  the  expected  peak  time.  Setting  r=0  in  (63)  and  (64)  we  get: 
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After  some  manipulations,  it  is  possible  to  verify  that: 
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Figs.  30  and  31  show  the  zero-Doppler  and  the  zero-delay  cuts  of  the  AAF  in  Fig.  29, 
respectively.  Recalling  that  the  definition  of  the  Rayleigh  resolution  is  the  peak-to-first  null 
distance,  by  inspection  of  Figs.  30  and  31,  it  is  apparent  that  the  delay  resolution  is  T  and  the 
Doppler  resolution  is  1/AT.  Considering  that  the  carrier  frequency  of  a  UMTS  base  station  if 
fc  =2100  MHz,  we  have  that  the  range  resolution  of  the  considered  waveform  is  of  39  m  and 
the  velocity  resolution  is  of  almost  0.71  m/s.  These  results  are  in  agreement  with  the  empirical 
results  showed  in  [Pet09],  [Gri03]  and  [Gri05],  where  the  authors  measured  and  analyzed  the 
AAF  of  a  UMTS  signal. 
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Figure  30  -  Zero-Doppler  cut  of  the  AAF  a  QPSK  signal  with  NT=0.1s,  T=0.26ps,  a=0.22. 


Figure  31  -  Zero-delay  cut  of  the  AAF  of  a  QPSK  signal  with  AT=0.1s,  7=0.26ps,  a= 0.22. 

In  monostatic  configuration,  estimation  of  the  time-delay  and  Doppler  shift  directly 
provides  information  on  target  range  and  velocity.  As  discussed  in  previous  sections,  this 
information  can  be  retrieved  also  in  a  bistatic  radar  configuration,  even  if  the  relations 
between  measured  delay  and  Doppler  frequency  and  target  distance  and  velocity,  respectively, 
are  not  linear  and  depend  on  the  geometry.  As  previously  shown,  referring  to  the  bistatic 
geometry  of  Figure  10,  the  delay  and  the  Doppler  shift  should  be  obtained  using  the  following 
geometry  dependent  non-linear  equations: 


-54- 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


(68) 


Rr  +  -\j  Rr  +  L  +  2Rk  L  si  n  0R 

c 


2—Vb  I -  + 


1 


Rr  +  L  sin  0R 


2  2^RR+l}  +  2RRLsm0R 


(69) 


Similar  relations  hold  for  the  actual  delay  xa  and  the  actual  Doppler  shift  va.  To  obtain  the 
Bistatic  Average  Ambiguity  Function  (BAAF)  in  the  range-velocity  plane  we  must  calculate 
x —tH-ta  and  v=vn-va  and  substitute  them  in  (61).  Therefore,  it  is  apparent  that  in  the  bistatic 
configuration,  the  BAAF  depends  on  the  transmitted  waveform  and  on  the  bistatic  geometry 
parameters,  i.e.  the  target  Direction  Of  Arrival  (DOA),  the  bistatic  baseline  length,  and  the 
distance  between  the  target  and  the  receiver. 

Figs.  32  and  33  show  the  cuts  of  the  BAAF  along  the  range  and  the  bistatic  velocity, 
respectively.  These  figures  have  been  obtained  assuming  V«=5  m/s,  Ra= 3  Km  and  L=5  Km  for 
three  different  values  of  the  receiver  look  angle  Or.  It  is  evident  that  in  all  the  cases,  the 
BAAF  reaches  its  maximum  when  the  hypothesized  range  and  velocity  equal  the  actual  range 
and  velocity  of  the  target.  Anyway,  the  shape  of  the  BAAF  strongly  depends  on  the  target 
angle  Or.  The  worst  case  is  Or  =— tt/2,  that  is,  when  the  target  is  on  the  baseline.  If  the  target  is 
between  the  transmitter  and  the  receiver,  the  BAAF  is  flat  and  the  range  and  velocity 
resolutions  are  completely  lost.  This  can  be  appreciated  by  realising  that  in  this  case  the  echo 
arrives  at  the  receiver  at  the  same  instant  as  the  direct  signal,  independent  of  the  target 
location,  and  the  Doppler  shift  of  a  target  crossing  the  bistatic  baseline  must  be  zero,  because 
the  transmitter-to-target  range  changes  in  an  equal  and  opposite  way  to  the  target-to-receiver 
range,  independent  of  the  magnitude  and  direction  of  the  target  velocity.  Note  also  that  for 
values  of  Or  far  from  -ji/2  the  shapes  of  the  BAAF  are  basically  the  same  and  their  behaviour 
is  similar  to  that  obtained  for  a  monostatic  system. 
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Figure  32  -  BAAF  as  a  function  of  range,  14=5  m/s,  74=3  Km  and  L=5  Km. 
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Figure  33  -  BAAF  as  a  function  of  the  bistatic  velocity,  Va=5  m/s,  7?«=3  Km,  L=5  Km. 


8.3  Monostatic  and  Bistatic  CRLB  of  the  UMTS  signal 


Similarly  to  previous  Section,  here  we  calculate  the  CRLB  exploiting  the  relationship 
between  the  AF  and  the  Fisher  Information  Matrix  (FIM)  of  the  delay  and  Doppler  parameters 
derived  in  [Van71]: 


Jc(t,v)  =  -2SNR- 


d2\zc{wf 

d2\zc(w)  f  1 

dr2 

drdv 

d2\zAwf 

d2\zc(w)\2 

dvdr 

dv2  \ 

(70) 
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In  (70)  we  use  the  subindex  c  to  highlight  that  this  relation  holds  for  deterministic  signals  and 
hence  it  involves  the  AF  for  a  fixed  sequence  of  symbols  c.  Therefore,  in  our  case,  the  FIM  in 
(70)  is  actually  a  conditional  FIM  (conditioned  to  a  given  sequence  of  symbols  c).  The 
modified  CRLB  (MCRLB)  immediately  follows  from  the  conditional  FIM,  since  the  MCRLB 
is  given  by  the  inverse  of  the  modified  FIM  (MFIM),  which  is  obtained  by  averaging  the 
conditional  FIM  of  (70)  with  respect  to  the  symbols  c  [Gin98],  [GinOO].  Hence,  by  averaging 
w.r.t.  c  the  conditional  FIM  in  (70)  and  then  taking  its  inverse,  we  get  the  MCRLBs  on  the 
estimation  of  target  time-delay  and  the  Doppler  shift  (i.e.  the  monostatic  MCRLBs).  It  is 
worth  remembering  that  to  calculate  the  classical  CRLB  on  time-delay  and  the  Doppler  shift 
we  should  average  out  the  dependence  of  the  random  symbols  c  in  the  joint  probability 
density  function  (pdf)  of  received  data  and  symbols.  This  is  infeasible  in  our  case.  Hence  we 
resort  to  calculation  of  the  MCRLB,  which  is  much  easier  to  obtain.  In  fact,  it  is  derived  from 
the  conditional  FIM  of  (70)  as  previously  outlined.  The  MCRLB  was  derived  in  [Gin98]  and 
its  relation  to  the  CRLB  investigated  e.g.  in  [Gin98],  [GinOO]  and  [Noa09],  and  it  was  found 
that  the  MCRLB  is  looser  than  the  CRLB.  The  MCRLB  was  applied  in  radar  e.g.  in  [Gin98b], 
[GinOOb],  and  [For  11]. 

To  evaluate  the  MFIM  and  hence  the  MCRLBs,  first  note  that  the  elements  of  the  MFIM 
can  be  expressed  as  [Van71]: 


J(gv) 


£{je(Gv)}  =  25AR 


4  7T2E[f2-f2]  2 xE{ft-ft\ 
2  7TE\Jt-Jt]  E{t2~~t2} 


(71) 


where: 


J  f2\Uc(f)\2df 

+°°  ~  I 

jfpAffdf  \. 

f  +°o 

|  J  t2  \uc(t)\2  dt 

+oo 

j*  t  |we  (t)|2  dt  L 


(72) 


(73) 


(74) 


(75) 
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(76) 


E{f} 


■  E  \  Im  < 


du  (?)  * .  .  .  . 

t—r - U  (t)dt 

c)t 


[  +oo  +00 

E[f~t)  =  E\  J  f\Uc(f)\2df  \t\uc(tfdt 


(77) 


The  expectation  E\ • }  is  with  respect  to  the  symbols  c,  Im{*}  is  the  imaginary  part  operator 
while  Uc(f)  is  the  Fourier  Transform  of  uc(t)  of  eq.  (50);  again  we  stress  that  we  used  the 
subindex  c  to  highlight  that  the  transmitted  signal  is  considered  for  a  fixed  sequence  of 
symbols  c.  Note  that  also  the  elements  of  the  MFIM  matrix  defined  in  (71)  directly  measure 
the  curvature  of  the  transmitted  signal  in  the  time-frequency  domain.  In  [TR5]  we  calculated 
the  elements  of  the  MFIM  in  (71)  using  the  definition  of  the  QPSK  signal  in  (71),  the 
obtained  results  show  that: 


J{t,v)  =  2SNR- 


-^—(x2  +3  arn2 
3  T2  v 

0 


24a2)  0 

T2(  1  N2-!) 

- - 1 - 

4  ^4#  3  J 


(78) 


The  MCRLBs  are  given  by  the  inverse  of  the  MFIM  of  (78),  being  J(r,v)  diagonal,  we  have: 


MCRLB(r)  =  -p — - - 

[J(^)]u 


MCRLB(v) 


1 


(79) 

(80) 


In  the  monostatic  configuration,  we  have  r  =2 R/c  and  v=2  Vfc/c,  therefore  it  is  easy  to  verify 
that  the  MCRLBs  of  the  estimated  parameters  are: 


MCRLB(R) 


_ 3T2c2 _ 

SSNR  [k2  +  3a2 7t2  -  24a2 )  ’ 


(81) 


MCRLB(V) 


c2 

T2 

(  i  a2-iV 

_l_ 

8  f2SNR 

4 

v4a  3  J 

(82) 


When  SNR=0dB,  NT=0A  s,  7=0.26  ps,  a=0.22  and  /c=2100MHz, 
the  values  of  the  root  MCRLB  (RMCRLB)  are  given  by:  RMCRLB(R)~15  m  and 
RMCRLB(V)  ~1.75  m/s.  Note  from  (79)  and  (81),  that  the  MCRLB  of  the  range  is  inversely 
proportional  to  [J(r,v)]i,i,  that  indicates  the  effect  of  the  signal  shape.  As  apparent  from  (71), 
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this  term  is  a  direct  measure  of  the  effective  bandwidth.  This  is  logical  because  increasing  the 
signal  bandwidth,  that  is,  the  value  of  a,  the  signal  has  a  faster  rise  time  giving  a  lower 
MCRLB  of  the  target  range.  On  the  other  hand,  from  (80)  and  (82)  we  have  that  the  bound  on 
the  velocity  estimation  accuracy  is  determined  by  [J(r,v)]2,2  that,  as  apparent  from  (71), 
measures  the  effective  pulse  length.  As  a  matter  of  fact,  the  larger  is  the  number  of  integrated 
pulses  N,  the  lower  is  the  MCRLB  on  the  target  velocity  estimation. 

As  done  in  previous  sections,  for  the  calculation  of  the  bistatic  MCRLB s  we  can  use  the  result 
on  the  monostatic  conditional  FIM  in  (70)  using  the  “chain  rule”  of  the  derivative.  In 
particular,  the  bistatic  conditional  FIM  is  given  by: 

•M*«V„)  =  PJ.(r.OPr,  («) 


where,  as  in  the  previous  sections, 


dr  dv 


dRR  c)Rr 
dr  dv 


3P.  SVB 


(84) 


Finally,  by  averaging  with  respect  to  c  we  get  the  bistatic  MFIM  as: 

J„(S„.Vs)  =  £{j,f(/;B.V's)}  =  P£{j1,(r,|/)}Pr=PJ(r.OPr  (85) 


The  relationship  in  eq.  (85)  shows  that  matrix  P  takes  into  account  only  the  effect  of  the 
bistatic  geometry  while  matrix  J(r,v)  takes  into  account  the  effect  of  the  transmitted  waveform 
and  the  SNR.  Note  that  also  the  SNR  depends  on  the  geometry,  in  fact  it  is  inversely 
proportional  to  the  energy  path  loss  factor  (Rr-Rt)2  due  to  propagation.  From  eq.  (85),  the 
bistatic  MCRLB s  of  the  receiver-to-target  range  and  the  bistatic  velocity  are  given  by: 


MCRLB  (Rr)  = 


[j»L 

[j.]u[j»L-[j»t 


(86) 


MCRLB(P„)  = 


(87) 


Note  that  the  monostatic  configuration  is  a  particular  case  of  the  bistatic  case,  and  it  is 
obtained  for  baseline  L= 0.  In  this  case  P  is  a  diagonal  matrix  whose  elements  are  [ P ]  u  =2/c 
and  [P]2,2=2-/c/c,  and  the  MCRLBs  of  range  and  velocity  are  the  same  as  in  eqs.  (81)-(82), 
respectively. 
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Figures  34-37  show  the  root  MCRLBs  (RMCRLBs)  of  range  and  velocity,  both  in  the 
monostatic  and  bistatic  case.  These  bounds  have  been  obtained  by  considering  the  QPSK 
signal  emitted  by  a  UMTS  base  station,  that  is  setting  T=  0.26  (.is,  a=0.22, /c=2100MHz  and 
NT=0.l  s.  All  these  Figs,  have  been  obtained  for  VB= 5  m/s  and  maintaining  the  SNR  at  OdB,  to 
highlight  only  the  effects  of  the  geometry.  In  particular,  Fig.  34  shows  the  RMCRLBs  as  a 
function  of  the  receiver-to-target  range  Rr,  for  two  values  of  the  receiver  look  angle, 
6r=- 0.4971  and  0K=n.  For  the  bistatic  system  the  baseline  has  been  set  to  5  km.  It  is  evident 
that,  for  all  the  parameter  values  we  investigated,  the  bistatic  RMCRLBs  are  always  higher 
than  the  monostatic  ones.  When  the  distance  from  receiver  to  the  target  increases,  the  bistatic 
system  behaves  more  and  more  as  the  monostatic  one. 

As  evident,  the  effects  of  geometry  are  prominent  where  the  target  approaches  the 
baseline,  that  is  in  the  case  of  0R  =-0.49n  when  Rr  <L.  In  this  case  the  RMCRLBs  tend  to  be 
very  high.  Similar  considerations  can  be  drawn  from  Fig.  35,  where  the  RMCRLBs  are 
plotted  as  a  function  of  the  baseline  L.  The  effects  of  the  geometry  are  more  apparent  in  Fig. 
36,  that  shows  the  RMCRLBs  as  a  function  of  the  receiver  look  angle  0R  ,  both  in  the  case  of 
Rr<L  and  Rr>L.  In  the  case  of  Rr<L.  when  the  target  approaches  the  baseline  (i.e.  0R  — >-7t/2), 
the  RMCRLBs  tend  to  infinity  because  the  resulting  delay  is  always  L/c  independent  of  the 
target  position  along  the  baseline,  and  the  radial  velocity  is  always  zero  independent  of  the 
target  speed  and  direction.  This  behaviour  is  in  agreement  with  the  shape  of  the  BAF  plotted 
in  Figs.  32-33. 

Fig.  37  shows  the  RMCRLBs  as  a  function  of  the  observation  time  NT.  As  expected  only 
the  RMCRLB  of  the  target  velocity  depends  on  this  parameter  and  decreases  as  the 
observation  time  increases.  Note  that  when  the  target  is  close  to  the  baseline,  the  RMCRLB  of 
the  velocity  is  very  high  also  for  increasing  values  of  the  observation  time.  On  the  other  hand, 
when  the  target  is  far  from  the  baseline  the  performance  of  the  monostatic  system  is  almost 
the  same  of  the  monostatic  one. 
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Figure  34  -  Root  MCRLBs  as  a  function  of  range  Rr,  NT=0. 1  s,  7=0.26  gs,  a=0.22, 
/c=2100MHz,  Va=5  m/s,  L=5Km,  SNR=0dB.  Monostatic  and  Bistatic  case. 
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Figure  35  -  Root  MCRLBs  as  a  function  of  the  baseline  L,  NT=0.l  s,  7=0.26  gs,  a=0.22, 
/c=2100MHz,  Va=5  m/s,  Rr= 3Km,  SNR=0dB.  Monostatic  and  Bistatic  case. 
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Figure  36  -  Root  MCRLBs  as  a  function  of  receiver  look  angle  0R,  NT-0. 1  s,  7=0.26  gs, 
a=0.22,  f(  -2 1 00MHz,  Ffl=5  m/s,  L=5Km,  SNR=0dB.  Monostatic  and  Bistatic  case. 
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Figure  37  -  Root  MCRLBs  as  a  function  of  the  observation  time  NT,  T- 0.26  gs,  a-0.22, 
fc= 2100MHz,  Vo=5  m/s,  Rr= 3Km,  L=5Km,  SNR=0dB.  Monostatic  and  Bistatic  case. 
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9  Sensor  selection  for  target  tracking 

This  Section  deals  with  the  problem  of  sensor  selection  of  the  optimum  channel  in  a 
multistatic  radar  system  along  the  trajectory  of  the  tracked  target. 

To  this  end  we  will  calculate  the  Posterior  Cramer-Rao  Lower  Bound  (PCRLB)  for  a  bistatic" 
radar  system.  The  PCRLB  is  a  useful  tool  which  bounds  the  error  variance  of  the  sequential 
target  state  estimates  produced  from  the  radar  measurements. 

The  Section  is  organized  as  follows.  Section  9.1  presents  the  mathematical  formulation  of  the 
problem  and  introduces  the  mathematical  model  of  the  target  state  and  measurements 
equations  for  bistatic  radar  systems.  Section  9.2  derives  the  theoretical  PCRLB  for  bistatic 
target  tracking  both  in  the  additive  process  noise  case  and  in  the  zero  process  case.  While 
Section  9.3  describes  how  to  exploit  the  PCRLB  for  sensor  selection  in  a  multistatic  radar 
system. 


9.1  Target  state  and  measurements  signal  model 

The  objective  of  target  tracking  is  to  estimate  the  target  trajectories  of  a  moving  target.  A 
point-like  target  motion  model  describes  the  evolution  of  the  target  state  with  respect  to  time. 
Let  assume  that  the  target  motion  equation  is  described  by  the  following  dynamic  target  state 
corrupted  by  additive  noise  model  [Tay79]: 

**+i  =fA-(x^aJ  +  n,  (88) 

where  x*,  a*  and  n*  are  the  target  state,  the  target  acceleration  and  the  process  noise  sequences 
respectively  while  f*  is  a  known  function  of  the  state  vector.  The  state  vector  at  time  index  k  is 
defined  as  x*=  [.*;*,  xk ,  y*,  yk  ]T  where  (xk,  yi)  indicate  the  position  of  the  target  while  ( xk ,  yk ) 
are  the  target  velocity  component  along  x  and  y  axes  of  the  Cartesian  coordinate  system.  The 
acceleration  sequence  is  defined  as  a/~[  xk ,  yk  \T,  where  xk  and  yk  are  the  component  of  the 

target  acceleration.  Here  this  sequence  is  assumed  known  and  is  used  as  a  control  input 
parameter  to  describe  a  particular  target  trajectory. 

The  sequence  models  the  process  noise  that  caters  for  any  mismodeling  effects  or 
unforeseen  disturbance  in  the  target  motion  model.  In  this  work  we  analyzed  the  linear  case, 


2 


This  bound  is  also  valid  for  monostatic  configuration  and  can  be  obtained  fixing  the  baseline  length  to  zero. 


-65  - 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


that  is  the  case  in  which  f*  is  a  linear  function  of  the  state  vector.  In  particular,  we  assumed 
that 

f*K.a*)  =  Fxt+Gat  (89) 


where 

"1  T  0  0 

0  10  0 

F  = 

0  0  1  T 
0  0  0  1 


(90) 


Gr 


r2/ 2  TOO 
0  0  r2/ 2  T 


(91) 


and  T  is  the  radar  sampling  interval. 

Note  that  f*  is  linear  and  its  Jacobian  is  equal  to  F. 

We  also  assume  that  the  random  noise  sequence  n .*  is  zero-mean  white  Gaussian  with 
covariance  matrix: 


T3/  3  T2/ 2  0  0 

r2/ 2  TOO 

0  0  T3/  3  r2/2 

0  0  T1/ 2  T 


(92) 


where  g  is  a  control  parameter  which  takes  into  account  the  power  of  the  process  noise. 

The  objective  of  target  tracking  is  to  recursively  estimate  the  target  state  from  a  set  of 
measurement  zk: 

Z*=MX*)  +  w*  (93) 


where  h*  is  a  known  function  of  the  state  vector  and  w/,  is  a  measurement  noise  sequence. 

For  the  bistatic  tracking  problem  analyzed  in  this  work,  the  available  measurements  at  time  k 
are  the  range  from  receiver  to  target  and  the  bistatic  velocity  .  In  particular,  the  function  h* 
does  not  depend  on  k  and  is  given  by  h(x*)=[r*,  Vk\r=[_hr(xk),  hv(xk)\r. 

To  give  explicit  expression  of  h(x^),  it  is  easy  to  verify  that: 

(94) 


3  Whit  respect  to  previous  sections  and  for  easy  of  notation,  here  the  symbol  r  denotes  the  range  from  receiver  to 
target,  rt  the  range  from  transmitter  to  target  while  v  denotes  the  bistatic  velocity. 
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(95) 


where 


K  (**)  =  L 


\K  +  ykyk 

\l*k  +yl 


h  = - xk  ~  xk  > 

L  +  rk+rtk 


yk  =  — — —  {yk  +rk)~  yk 


L+rk  +  rh 


rtk  =  ^lxk+(yk~Ly 


(96) 

(97) 

(98) 


As  apparent  from  eq.  (93),  the  bistatic  measurements  are  affected  by  additive  noise  w*, 
independent  of  the  process  noise  n*.  The  additive  noise  w*  is  assumed  to  be  Gaussian 
distributed  with  zero  mean  vector  and  covariance  matrix  R/  . 

As  previously  shown,  in  bistatic  radar  systems  the  accuracy  of  the  estimate  of  range  and  the 
bistatic  velocity  heavily  depends  on  the  transmitted  waveform  and  on  the  geometry  of  the 
scenario.  In  previous  sections  it  was  shown  that  the  Fisher  Information  Matrix  (FIM)  of  the 
radar  measurements  can  be  expressed  as  Jf  =PA  JMP/,  where  JM  is  the  monostatic  FIM  of 
time-delay  and  Doppler  shift  and  depends  only  on  the  Ambiguity  Function  (AF)  of  the 
transmitted  waveform.  On  the  other  hand,  P/  depends  on  the  geometry,  which  changes  with 
time  k  while  the  target  is  moving  along  its  trajectory. 


9.2  PCRLB  for  bistatic  radar  systems 

Let  indicate  with  J \  the  filtering  information  matrix  at  time  k.  Its  inverse  is  the  PCRLB 
that  bounds  the  error  variance  of  the  target  state  estimate  at  time  k,  that  is 

J*1  ^  E{(x»  ~xkf}  (99) 

where  xklk  is  an  unbiased  estimator  of  the  state  vector  based  on  all  the  available  measurements 
up  to  time  k. 

In  [Tic98]  the  authors  provided  an  elegant  method  of  recursively  computing  the  filtering 
information  matrix  J*.  In  particular,  when  the  additive  noise  sequences  n/t  and  Wk  in  eqs.  (89) 
and  (93)  are  Gaussian  and  mutually  independent  and  their  covariance  matrix  are  invertible, 
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the  information  matrix  can  be  computed  recursively  using  the  following  equation 
[Tic98],[Tay79]: 

J»+,=(Q+rji-,Fr)"'+E{HLR-;iHltl}  (100) 

where  H^+i  is  the  Jacobian  of  h(x*+i)  evaluated  at  the  true  state  x*+i;  that  is: 

dh  dh.  dh.  dh 
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&-I-1 


dx,_,  dx, 


xk+ 1 
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dh 
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'  k+l 


(101) 


The  recursive  equation  in  (100)  can  be  viewed  as  the  sum  of  two  terms,  the  first  one 
(Q+FJi  F7)"1  represents  the  a  priori  information  given  by  the  previous  target  state,  while  the 
second  term  E{  }  represents  the  information  gained  by  the  radar  measurements. 


The  expressions  of  the  elements  of  H*+i  are  given  by 
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where  dk+l  =  yjxl+i  +  vA2+1  .  The  derivatives  that  appear  in  (104)  and  (105)  can  be 

straightforwardly  calculated  and  are  not  reported  here  for  lack  of  space.  From  these  equations 
it  is  clear  that  along  the  target  trajectory,  the  FIM  J*  depends  on  sensor  accuracy  (through  R/., 
which  is  itself  dependent  on  the  target  trajectory  and  on  the  transmitted  waveform),  the 


-68- 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


sampling  interval  T  and  the  baseline  length  L.  In  particular,  as  shown  in  the  previous  section, 
when  the  target  is  crossing  the  baseline,  resolution  is  totally  lost  and  therefore  the  errors  of  the 
measurements  tends  to  infinity.  Matrix  1  tends  to  zero  and  hence  also  the  second  term  in 
(100).  In  this  case  there  is  no  information  gain  collected  by  the  target  measurements. 
Regarding  eq.  (100),  the  most  critical  problem  in  calculation  appear  to  be  the  expectation 
operator  which  is  taken  with  respect  to  the  state  vector  x*.  This  calculation  can  be 
approximated  through  Monte  Carlo  simulation  by  generating  a  set  of  state  vector  realizations 
and  then  computing  the  average  over  this  set.  However,  depending  on  the  analyzed  problem, 
this  is  rarely  done  in  practice.  Instead,  often  purely  deterministic  trajectories  are  considered. 
When  the  target  motion  equation  is  considered  in  absence  of  process  noise,  the  PCRLB  refers 
to  the  zero  process  noise  case.  Clearly,  in  this  case  the  obtained  PCRLB  is  a  more  optimistic 
lower  bound.  In  absence  of  process  noise,  the  evolution  of  the  target  state  is  purely 
deterministic  and  hence  the  expectation  operator  in  eq.  (100)  can  be  dropped  out.  Moreover, 
the  covariance  matrix  Q  becomes  zero  and  therefore  it  follows  that  [Tay79]: 

J<«=[F_,]rJ,F-'+H[tlR1-;iH„  (107) 

Clearly,  also  in  the  zero  process  noise  case,  the  filtering  information  matrix  is  the  sum  of  two 
terms,  one  takes  into  account  the  a  priori  information  while  the  other  the  measurement 
information.  For  L= 0  the  bistatic  PCRLB  s  coincide  with  the  monostatic  PCRLB  s. 


9.3  Application  of  PCRLBs  to  sensor  selection 

In  this  Section  we  show  how  the  PCRLBs  can  be  used  to  sensor  selection  in  a  multistatic 
radar  system.  In  particular  we  consider  a  system  in  the  context  of  maritime  border  control 
scenario,  composed  by  a  Passive  Coherent  Location  (PCL)  receiver  that  exploits  the  signal 
emitted  by  two  non  co-operative  transmitters  of  opportunity:  a  UMTS  base  station  and  a  FM 
radio  station.  These  kind  of  systems  have  some  significant  attractions  because  it  allows  the 
use  of  parts  of  the  RF  spectrum  (VHF/UHF)  and,  since  broadcast  transmissions  at  these 
frequencies  can  have  substantial  transmit  powers,  they  can  have  excellent  coverage. 
Furthermore,  advantage  should  be  drawn  utilising  the  UMTS  sources  of  illumination.  In  fact, 
the  FM  system  could  detect  targets  at  relatively  far  range  and  use  its  fine  Doppler  resolution 
to  build  a  crude  track.  Once  the  target  appears  within  the  coverage  area  of  the  UMTS  system, 
this  can  be  exploited  for  its  superior  range  resolution  properties.  The  aim  of  this  section  is  to 
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evaluate  the  performance  of  each  bistatic  channel  in  estimating  the  target  trajectory  of  a  target 
approaching  Leghorn  harbour  and  also  to  propose  an  algorithm  for  optimally  selecting  the 
transmitter  of  opportunity.  This  algorithm  can  provide  a  significant  aid  to  harbour  protection 
and  can  reduce  the  computational  load  of  surveillance  operations. 

The  analyzed  multistatic  scenario  is  shown  in  Figure  38.  The  receiver  is  placed  in  “via 
Michelangelo”,  between  “darsena  Petroli”  and  “darsena  Pisa”  of  Leghorn  harbour.  Its  antenna 
has  a  gain  G=10dB  and  an  Half  Power  Beam  Width  (HPBW)  of  3°.  It  exploits  the  signals 
emitted  by  two  non  co-operative  transmitters.  The  first  one  is  a  UMTS  Base  Station  located  in 
“Piazzale  Marmi”,  453m  away  the  receiver  in  the  South-East  direction.  The  second 
transmitter  is  an  FM  commercial  radio  station  located  in  “Monte  Serra”,  36Km  away  the 
receiver  in  the  North-East  direction. 
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Figure  38  -  Analyzed  multistatic  radar  System. 
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Figure  39  shows  the  Root  CRLB  (RCRLBs)  of  target  range  and  bistatic  velocity,  that  is  the 
root  of  the  diagonal  elements  of  Ra.  The  RCRLBs  of  the  UMTS  channel  are  shown  at  the  top 
of  the  figure,  while  the  RCRLBs  of  the  FM  channel  are  shown  at  the  bottom.  To  calculate 
these  maps  we  supposed  that  the  target  was  moving  horizontally  from  left  to  right  with  a 
constant  velocity  of  3m/sec.  From  the  results  it  is  clear  that  the  performance  in  measuring  the 
target  range  and  the  bistatic  velocity  strongly  depends  both  on  the  geometry  and  the 
transmitted  waveform.  It  is  also  interesting  to  note  that  the  performance  in  estimating  the 
target  velocity  is  better  in  the  FM  channel  than  in  the  UMTS  one,  due  to  the  longer  time  of 
observation.  On  the  other  hand,  the  range  resolution  is  better  in  the  UMTS  channel  than  in  the 
FM  one.  Note  also  that  the  FM  channel  operates  with  very  high  transmitted  power  giving  a 
better  coverage  and  appreciable  resolution  also  in  the  far  range.  Using  the  information 
provided  in  Figure  39  combined  with  the  results  obtained  in  previous  sections,  it  is  possible  to 
select,  for  each  point  of  the  analyzed  area,  the  channel  having  the  best  performance.  For  each 
location  of  the  surveillance  area  it  is  possible  to  calculate  the  measurement  information  matrix 
jf+1  =  H^+1R^,Ha+1  and  then  to  provide  a  function  of  this  matrix  that  evaluates  the  channel 

that  gives  the  higher  measurement  information.  As  an  example,  let  consider  the  case  of  a 
naval  target  approaching  the  harbour  with  constant  velocity  equal  to  3  m/sec  (see  the  route 
Palermo-Livorno  (Leghorn)  of  Fig.  38).  When  the  target  trajectory  is  deterministic  the  bistatic 
PCRLBs  of  each  channel  are  obtained  by  inverting  matrix  Ji+1  of  eq.  (107).  Note  from  (107) 

that  the  target  dependent  quantity  det(  R l:  1 )  can  be  considered  as  an  index  of  the  measurement 
information  and  hence  can  be  used  to  select  the  channel  with  the  best  performance. 

We  investigated  also  the  performance  that  can  be  obtained  using  different  indices  of  the 
measurement  information,  such  as  the  trace,  the  determinant  and  the  norm  of  the  matrix 
H^+jR^H^j  .  The  obtained  performance  are  similar,  anyway  the  best  performance  are 

achieved  using  the  index  det(  R^1 ).  This  is  justified  considering  that  for  the  particular  case  we 

analyzed,  the  parameters  that  describe  the  uncertainty  on  the  state  and  the  measurement 
vectors  are  only  the  probability  of  detection  and  measurement  error  covariance  matrix  and 
hence  the  index  of  the  measurement  information  must  be  a  function  of  these  parameters. 
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RCRLB  Range  -  Channel  1 


RCRLB  Velocity  -  Channel  1 


Figure  39  -  Root-CRLB  of  range  (left  hand  side)  and  bistatic  velocity  (right  hand  side) 
for  the  UMTS  channel  (top)  and  FM  channel  (bottom). 


-72- 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


[dBm/sec]  [dBm/sec] 


Best  Channel 


Figure  40  shows  the  channel  which  has  the  higher  information  for  each  point  of  the  analyzed 
area.  The  scale  of  this  figure  is  quantized  into  2  levels:  level  1  refers  to  the  UMTS  channel 
while  level  2  to  the  FM  channel.  In  Figure  40  we  also  plotted  a  red  line  that  represents  the 
trajectory  of  the  ship  approaching  the  harbour.  Being  the  target  trajectory  deterministic,  the 
PCRLBs  of  the  receiver  that  dynamically  select  the  best  channel  are  given  by: 

J,«=[F"]rJlF-1+J“  (108) 


where: 

|[H™]r[RS,]"H®,  ifdet([R™  ]"j>det([RS]"j 
[H®,  I  [»£’,  ]"  H£>  if  del  ([R™,  J  j  <  del  ([R™]- j 


(109) 


that  is,  at  each  step  k+ 1  the  measurement  information  term  J™  is  given  by  selecting  the 

channel  with  the  higher  measurement  information.  In  (109)  the  upper  index  (1)  refers  to  the 
UMTS  channel  while  the  index  (2)  to  the  FM  channel.  The  PCRLBs  of  the  receiver  that  select 
the  best  transmitter  are  a  shown  in  Figure  41,  where  we  also  plotted  the  PCRLBs  of  each 
single  channel.  From  the  results  it  is  evident  that  there  is  a  substantial  gain  with  respect  to 
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each  bistatic  channel  and  the  obtained  performance  is  equal  or  one  order  of  magnitude  better 
than  the  performance  of  the  channel  with  the  lowest  PCRLB.  As  apparent,  initially  the 
performances  are  the  same  as  the  FM  channel.  This  is  due  to  the  fact  that  in  the  far  range  the 
FM  channel  has  the  best  performance  thanks  to  the  higher  SNR.  When  the  target  approaches 
the  harbour,  the  UMTS  channel  has  a  better  resolution  and,  exploiting  this  channel,  the 
proposed  receiver  is  able  to  increase  the  performance  in  estimating  the  target  trajectory. 
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Figure  41  -  Root-PCRLBs.  UMTS  channel,  FM  channel  and  dynamic  selection  channel. 
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Conclusions 

The  considerable  flexibility  of  modem  radar  systems  calls  for  new  methods  of  designing 
and  scheduling  the  waveforms  to  optimize  the  radar  performance  and  to  compensate  for  the 
environment  and  multi-radar  interference.  To  this  end  the  first  part  of  our  research  activity 
focused  on  the  design  and  analysis  of  a  set  of  waveforms  in  the  contest  of  netted  radar 
systems  that  guarantee  good  target  detection  and  parameter  estimation  in  different  scenarios 
and  allow  an  optimal  access  to  the  same  transmit  channel.  In  particular  we  focused  on  the 
analysis  of  the  Frequency  Coded  Waveforms  where  the  frequency  hop  patterns  were 
constructed  upon  an  extension  of  the  theory  of  quadratic  congruences. 

These  signals  are  characterized  by  an  auto-ambiguity  function  that  exhibits  a  narrow  thumb 
tack  shape  with  low  sidelobes  and,  in  contrast,  with  small  cross-ambiguity  functions.  In  our 
work  we  also  investigated  the  impact  of  an  interfering  signal  emitted  by  a  non  cooperative 
radar  on  the  performance  of  the  reference  radar  in  the  estimation  of  the  target  direction  of 
arrival  (DO A).  Our  research  also  focused  on  the  design  and  the  analysis  of  adaptive  waveform 
diversity  technique  to  mitigate  the  impact  of  the  interfering  radar  on  the  performance  of  the 
reference  radar.  To  this  end,  we  presented  four  DOA  estimators,  the  Modified  Pseudo- 
Monopulse  (MPM),  the  Modified  Maximum  Likelihood  (MML),  the  Sub-sequent  Pseudo- 
Monopulse  (SPM)  and  the  Sub-sequent  Maximum  Likelihood  (SML). 

Our  research  activity  on  waveform  diversity  and  design  also  concerned  on  Multistatic  active 
radar  systems,  that  is  systems  that  exploit  multiple  transmitter  and  receiver.  Such  systems 
differ  from  typical  modem  active  radar  through  consisting  of  several  different  monostatic  and 
bistatic  channels  of  observation.  Due  to  this  spatial  diversity,  these  systems  present  challenges 
in  managing  their  operation  as  well  as  in  usefully  combining  the  multiple  sources  of 
information  on  a  particular  area  of  surveillance. 

The  information  gain  obtained  through  this  spatial  diversity  combined  with  some  level  of 
Knowledge-based  data  processing  can  give  rise  to  a  number  of  advantages  over  both  the 
individual  monostatic  and  bistatic  cases  in  typical  radar  functions,  such  as  detection, 
parameter  estimation,  tracking  and  identification.  As  known,  the  performance  of  each  bistatic 
channel  heavily  depends  upon  the  geometry  of  the  scenario,  that  is  the  position  of  the  target 
with  respect  to  each  receiver  and  transmitter,  and  the  transmitted  waveform.  Therefore,  the 
main  goal  of  our  research  was  the  develop  of  mathematical  tool  that  can  be  used  to  evaluate 
the  best  achievable  performance  of  monostatic  and  bistatic  systems.  To  this  end  we  calculated 
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the  Ambiguity  Function  (AF)  and  the  Cramer- Rao  Lower  Bounds  (CRLB)  of  several 
transmitted  waveforms  for  both  the  monostatic  and  bistatic  configurations. 

Our  research  also  focused  on  the  analysis  of  passive  systems,  in  particular  we  analyzed  the 
performance  of  PCL  (Passive  Coherent  Location)  systems,  that  is  passive  bistatic  radar 
systems  where  the  transmitter  of  opportunity  is  from  a  non-radar  transmission,  such  as 
broadcast,  communications  or  radionavigation  signal.  These  systems  have  the  advantage  that 
the  receivers  do  not  need  any  transmitter  hardware  of  their  own,  and  are  completely  passive, 
and  hence  undetectable.  PCL  multistatic  systems  have  some  significant  attractions  because 
they  can  allow  the  use  of  parts  of  the  RF  spectrum  (VHF  and  UHF)  that  are  not  usually 
available  for  radar  operation,  and  which  may  offer  a  counterstealth  advantage,  since  stealth 
treatments  designed  for  microwave  radar  frequencies  may  be  less  effective  at  VHF  and  UHF. 
Moreover,  broadcast  transmissions  at  these  frequencies  can  have  substantial  transmit  powers 
and  the  transmitters  are  usually  sited  to  give  excellent  coverage.  In  particular,  in  our  research 
we  derived  the  AF  and  the  CRLBs  on  the  estimation  accuracy  of  the  target  range  and  velocity 
for  a  QPSK  signal  where  the  pulses  are  shaped  with  a  Root  Raised  Cosine  filter.  These  results 
provide  a  mathematical  tool  to  evaluate  the  performance  of  a  Passive  Coherent  Location 
(PCL)  system  where  the  receiver  exploits  the  downlink  signal  emitted  by  non  co-operative 
UMTS  base  station.  The  actual  growing  coverage  of  UMTS  signals  on  the  international 
territory  makes  a  multistatic  radar  configuration  feasible.  All  the  obtained  results  can  be 
extended  in  a  straightforward  manner  to  any  Pulse  Amplitude  Modulation  (PAM)  signal,  in 
order  to  analyze  the  performance  of  a  PCL  system  that  exploits  the  signal  emitted  by  another 
digital  non  co-operative  transmitter  of  opportunity. 

Our  research  activity  also  concerned  on  the  design  of  an  algorithm  for  optimally  selecting  the 
best  channel  (tranmitter-receiver  pair)  for  the  tracking  and  the  parameter  estimation  of  a  radar 
target  moving  along  a  certain  trajectory  in  a  multistatic  scenario.  To  this  end  we  exploited  the 
information  gained  through  the  calculation  of  the  CRLB.  The  idea  was  to  select  the  base 
station  that  gives  the  lowest  CRLB  for  the  target  velocity  or  range.  We  analyzed  as  well  the 
problem  of  the  optimum  sensor  selection  along  the  trajectory  of  a  tracked  target  in  a 
multistatic  radar  systems.  To  this  end  we  evaluated  the  Posterior  Cramer- Rao  Lower  Bound 
(PCRLB)  of  the  sequential  target  state  estimation  and  we  defined  an  algorithm  that  exploits 
this  mathematical  tool  to  select  the  best  channels  for  the  tracking  of  the  target.  Our  research 
also  focused  on  exploiting  the  CRLB  to  compute  the  rules  for  selecting  the  best  weighting 
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coefficients  for  fusing  the  signals  from  the  receivers  in  a  multistatic  system,  in  order  to 
improve  the  detection  of  the  target  and  the  estimation  accuracy  of  its  kinematic  parameters. 

It  is  important  to  consider  that  a  system  that  exploits  the  CRLB  to  compute  the  rules  for 
improving  its  own  performance,  is  a  system  that  has  “cognition”  and  changes  its  mode  of 
operation  according  to  the  way  in  which  it  perceives  its  operating  environments. 

Future  research  activity  on  waveform  diversity  and  design  can  be  focus  on  the  design  of 
Cognitive  radar  system.  The  increasing  needs  for  frequency  allocations  to  radiolocation 
services  is  in  contrast  with  the  growth  of  activities  in  the  area  of  civil  communications,  the 
emergence  of  new  technologies  and  new  services  that  involve  a  strong  demand  for  spectrum 
allocation  inducing  a  very  strong  pressure  upon  the  frequency  channels  currently  allocated  to 
radars.  For  this  reason,  the  availability  of  frequency  spectrum  for  multifunction  radar  system 
has  been  severely  compromised,  available  frequency  bands  have  been  diminished  and 
increasingly  dispersed. 

This  unique  issue  of  spectrum  crowding  and  steadily  increasing  radar  requirements  cannot  be 
addressed  by  traditional  modes  of  operation.  Multifunction  radar  systems  require  the  ability  to 
anticipate  the  behavior  of  radiators  in  the  operational  environment.  This  in  turn  leads  to  the 
need  for  critical  and  new  found  methodologies  based  upon  cognition  as  an  enabling 
technology. 

There  is  a  need  to  generate  situational  awareness  and  directly  support  the  delivery  of  effects  - 
then  intelligent  aggregation  of  sensor  data,  information  sources,  and  knowledge  in  real-time  is 
required.  Achieving  the  effects  of  cognition  is  essential  to  this  goal.  In  light  of  this  spectrum 
encroachment,  intelligent  spectrum  utilization  coupled  with  agile  waveform  selection  and 
signal  processing  in  multifunction  radars  has  become  a  necessity. 

The  concept  of  Cognitive  Radar  represents  a  way  forward  to  address  these  issues  from  a 
sensor  and  information  perspective.  The  Cognitive  Radar  will  achieve  the  effects  of  cognition 
in  the  context  of  radar  operation  using  perception,  memory  and  judgment.  The  Cognitive 
Radar  will  have  the  equivalent  effect  of  cognition,  i.e.,  ability  to  “reason,”  without  developing 
“consciousness”  through  the  effective  use  of  radar  signal/data/image/track  processing; 
waveform  generation,  timing  and  control;  platform  placement;  and  feedback  for  further 
adaptation  and  enhanced  control. 
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